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ABSTRACT 


ANALYSIS AND TEST OF A BREADBOARD CRYOGENIC 
HYDROGEN/FREON HEAT EXCHANGER 

by 

L. F. Desjardins and J. Hooper 
April 1973 


Ibis report describes the system studies required to verify a tube - in - 
tube cryogenic heat exchanger as optimum for the Space Shuttle mission. Design 
of the optimum configuration, which could be fabricated from commercially 
available hardware, is discussed. Finally ; testing of the proposed configuration 
with supercritical hydrogen and Freon 21 is discussed and results are compared 
with thermal and dynamic analysis. 
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FOREWORD 


This report was prepared by the Hamilton Standard Division of the 
United Aircraft Corporation for the National Aeronautics and Space Adminis- 
tration's Manned Spacecraft Center in accordance with Contract NAS 9-12725. 
This report covers the period from May 12, 1972 through January 24, 1973. 
During this period the tube - in - tube was selected as the optimum cryogenic 
heat exchanger configuration for the Space Shuttle. A breadboard unit was 
designed, fabricated and tested to correlate analysis with actual performance. 

Personnel responsible for the conduct of this program were Mr. F. H. 
Greenwood, Program Manager, Mr. L. F. Desjardins, Engineering Project Manager , 
and Messrs. J. Hooper and R. Balinskas, Analytical Engineers. Testing was 
conducted at Wyle Laboratories at Norco, California. Mr. F. Collier was 
Technical Monitor and Mr. R. J. Gillen overall program supervisor for the NASA 
Manned Spacecraft Center. 
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SUMMARY 


Initial studies indicate that the Space Shuttle system heat rejection 
requirements may be optimized by the use of a cryogenic heat exchanger. The 
eventual viability of this heat rejection method is dependent upon the 
technology development of a heat exchanger concept for use with cryogenic 
hydrogen and Freon 21. A simple tube - in -tube concept, selected for the 
Space Shuttle Orbiter under contract NAS 9-12208, showed promise of providing 
the required features at a very low weight penalty and providing basic 
conceptual simplicity. This concept had been modeled mathematically and its 
feasibility demonstrated by tests conducted with cryogenic nitrogen and water. 

The objective of thisphase II program was the correlation of analytically 
derived data with test data from a heat exchanger utilizing the Shuttle fluids 
of supercritical hydrogen and Freon 21. As can be seen in the Work Breakdown 
Structure (WBS) , figure 1 , the program was divided into 5 major tasks , which 
were sequential except for the generalized task of Management and Documentation 
which occurred throughout the program. At the outset, the latest available 
Shuttle data was examined and found similar to that on which the previous 
phase was based. Thus the same rationale supporting a tube - in - tube cryo- 
genic heat exchanger was still applicable, as well as the proposed control 
system. Next, the heat exchanger was optimally sized using commercially 
available tubing. Aluminum was selected because of its superior physical 
properties (1) at cryogenic temperatures , high conductivity and minimum weight 
The design selected, namely concentric 1" O.D. and. 3/4" O.D. concentric tubes , 
15 feet long in a coil, approximately 1 foot diameter, was included in the 
math model and computer program previously prepared. System dynamic studies 
showed the proposed cryogenic heat exchanger control system both accurate 
and stable. 

The proposed heat exchanger was defined in drawing SVSK 86021 and a 
vendor was selected for its manufacture after receiving NASA design approval . 
Controls of the proposed cryogenic heat exchanger control loop were not 
definitized in drawings since their procurement or test were not a part of 
the program. Considerable difficulty was encountered in manufacture of the 
heat exchanger in that the maximum length of commercially available tubing 
was 12 feet requiring a butt welded joint which fractured during coiling. 

In addition severe rippling of the outer tube was experienced. To correct 
these problems the wall thickness of the outer tube was increased from .028" 
to .035" and the tubing length was limited to 12 feet. Future prototype 
procurement should provide time and binding for a mill run of the proper 
length tubing. 

Wyle Laboratories in Norco, California was selected to conduct the 
heat exchanger tests after soliciting bids from other test facilities and 
examining Hamilton Standard's test capabilities. Wyle Was particularly 

(l)ALCOA Bulletin "ALCOA Aluminum The Cryogenic Metal." 
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advantageous in that they had a complete functional cryogenic hydrogen 
facility and took no deviations to the test plan. Steady state thermal 
performance tests were conducted at 5000, 20,000, 35,000, and 50,000 Btu/hr. 
Transient tests were conducted by increasing Freon inlet temperature from 
minimum to maximum in step and ramp changes and starting and stopping 
hydrogen flow at 5000 and 50,000 Btu/hr heat loads. Good correlation with 
analytically derived thermal and dynamic performance was obtained ; however, 
a progressively more severe hydrogen instability was noted with increasing 
heat load. This was manifest as a one hertz fluctuation in hydrogen inlet and 
outlet pressure and flow. Amplitude of the pressure fluctuations at the 
50,000 Btu/hr case was + 15-20 psi. Investigation confirmed that the hydrogen 
pressure/flow fluctuations were not caused by the test rig. 


Literature searches revealed a considerable number of similar instabil- 
ities but the mechanism of instability has been postulated only. Some show 
it as a function of low cryogenic flow rates and high heat loads and others 
as a function of cryogenic inlet to outlet density ratios. Clearly, however, 
the mechanism of instability must be understood better before using a system 
with the inherent instability or attempting to redesign a heat exchanger for 
stability. 
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INTRODUCTION 


This program was conducted in accordance with NASA Contract NAS 9-12725 
to verify the selection of a tube - in - tube configuration as the optimum 
cryogenic heat exchanger for the Space Shuttle program, and design, manufac- 
ture and test a breadboard unit. This program is a follow-on to NASA Contract 
NAS 9-12208 which authorized conduct of an in depth study to develop a cryogenic 
heat exchanger concept meeting Space Shuttle requirements of simplicity, light- 
weight, reusability, low maintenance, long life and low cost. This program 
recommended a tube - in - tube counterflow configuration on the basis of lower 
cost and longer life than the other primary contenders of tube - in - shell and 
plate and fin. System weight difference among these three primary contenders 
was only 3.6 of 121.4 pounds hence the small weight advantage of the plate 
and fin configuration could not counteract the cost and life advantages of the 
tube - in - tube confgiruation. 

The major tasks of this program consist of the following: 

Design 'Studies - These studies reviewed established results with latest 
Shuttle conditions and provide further math model development and system 
dynamic evaluation. 


Breadboard Unit Design - This task established a tube - in - tube config- 
uration which provided minimum system weight consistent with commercially 
available tubing and fittings. 

Breadboard Uttit Fabrication - A prototype heat exchanger was fabricated. 

Breadboard Tests - Steady state and transient tests were conducted within 
a range of 5000 to 50,000 Btu/hr to afford correlation with the analytical 
approach for dynamics and thermal performance. 
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CONCLUSIONS 


• The tube - in - tube heat exchanger is confirmed the optimum configuration 
for Shuttle cryogenic application. 

• Math model and computer studies have shown that the previously selected 
thermal control system provides the required system stability and 
accuracy. 

• Test data verifies thermal and dynamic analysis of the heat exchanger 
math model. 

• . The cryogenic heat exchanger as tested displays a hydrogen pressure 

and flow oscillation of approximately one hertz which increases in 
magnitude with increasing heat loads. 

• The thermal control system proposed for use with the cryogenic heat 
exchanger meets thermal stability and accuracy requirements despite 
hydrogen oscillations of the selected heat exchanger. 

• The mechanism of hydrogen oscillations is neither completely under- 
stood nor predictable. 
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REOOMM1ATIONS 


• Further studies and tests should be conducted to: 

Understand the mechanism of cryogenic hydrogen oscillations to 
a degree permitting stable heat exchanger design. 

Understand the characteristics and implications of such 
oscillations in a Shuttle cryogenic fluid loop. 

• Mill runs should be procured in the future to permit a heat exchanger 
length resulting in optimum system weight. 

• For the tubing sizes selected in this study coil size should be a 
minimum of 12 inches I.D. to prevent rippling. 
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DISCUSSION 


The discussion is divided into four major sections following the 
program work breakdown structure. These sections are Design Studies, Unit 
Design, Unit Assembly, and Unit Test and Delivery. 


DESIGN STUDIES 


A review of the cryogenic heat exchanger requirements and system penalties 
was couple ted. Contact with North American Rockwell showed that heat loads and 
cryogenic system penalties had not changed. Thus the tube - in - tube con- 
figuration determined best under NASA/MSC contract NAS 9-12208 continues as 
the selected configuration. An optimization of the tube - in - tube heat 
exchanger size and weight using commercially available aluminum tubing was 
conducted. This optimization considered a maximum heat load of 50,000 Btu/hr, 
and normal heat loads, zero load operation, hydrogen utilization, hydrogen 
supply pressure variations, degree of hydrogen super-heat and coolant freez- 
ing prevention. In addition, nominal loads were assumed for gross structural 
investigation as specified in figure 3. 

The cryogenic heat exchanger was optimized based on the following 
condition: 


T of Frecn 21 in 

= 

620°R 

T of Freon 21 out 

s 

500 °R . 

T of H 2 in 

= 

66 °R 

P of H 2 in 

= 

300 psi 

Q 

“ 

50,000 Btu/hr 

W of Freon 21 

= 

1602.5 lbs/hr 

Usage Time 

= 

1 hour 


This condition is based on North American Rockwell data. Power 
penalty to account for the Freon 21 pressure drop was calculated as 2.39 
lbs/psi based on the North American Rockwell system data. 

The pipe outside diameters that were used in the optimization were 
3/8" to 1" in commercially available 1/8" increments between these limits. 
Thicknesses investigated were .020, .028, and 035 inches. As aluminum is 
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compatible with both hydrogen and Freon 21 and has excellent mechanical 
properties at cryogenic temperatures , it was selected as the heat exchanger 
material. Due to the long tube length, conduction along the walls is 
insignificant. 

Figure 2 compares the following combinations of pipe and wall thiclaiess: 

1. 1/2", .028; 3/4", .028" 

2. 3/4", .028; 1", .028" 

3. 5/8", .028; 7/8", .028" 

4. 3/8", .020; 5/8", .020" 

5. 3/4", .035; 1", .035. 

Total equivalent weight consists of: 

a) Heat exchanger "wet weight". 

b) Hydrogen expended for one hour. 

c) Freon punping power penalty. 

d) Tankage penalty was not included. 

The difference in total equivalent weight for the configurations compared 
was less than two pounds at their optimum lengths. Configuration 2 was chosen 
with a length of 15 feet as it has the lowest total equivalent weight. 

All of the configurations examined had the minimum wall temperatures above 
410° R with the exception of configuration 4. As configuration 2 was selected 
there is no danger of freezing the Freon 21. Freon 21 freezes at 249°R. 

Variations in H 2 supply pressures do not affect the heat exchanger size 
and were not considered in heat exchanger sizing. 

For comparison, a plate fin heat exchanger was sized to determine the 
weight saving that could be achieved. A maximum weight saving of 3.6 pounds 
was calculated using a plate fin heat exchanger. It was concluded that this 
weight saving did not justify the increased cost and development required for 
a plate fin heat exchanger. 

The stresses of the coiled design were investigated for: 

a) Thermal stresses. 

b) Vibrational stresses. 

c) Shock. 
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CRYOGENIC DOUBLE PIPE HEAT EXCHANGERS 
FIGURE 2 
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The configuration, mounting positions, and requirements are shown in 
figure 3. The conclusion was reached that there are no structural problems 
in the design. 


Controls Simulation and Math Model Update 

The analytical model for the dynamics of a cryogenic tube - in - tube 
heat exchanger, which had previously been computerized, was updated to include 
the results of the optimization studies and was expanded to include the 
dynamics of the Shuttle Freon loop. 

The Freon loop parameters for this investigation were selected on the 
basis of recent information obtained from North American Rockwell Space 
Division as well as in-house studies. Figure 4 shows the design heat loads 
selected for the active components and loop temperatures . Figure 5 shows 
line lengths and the weight of all components in the loop. The dynamic 
parameters for the model are summarized in Table I. 

The computerized model includes all elements in the radiator bypass loop, 
but not the radiator, hydraulic heat exchanger and associated lines. The 
cabin water loop line from the interface heat exchanger to the cabin condensing 
heat exchanger has also been included in order to obtain the effect of Freon 
loop temperature fluctuations on the temperature of cooling water entering 
the condensor. 


System Dynamics Investigation 

A mathematical model of a tube - in - tube cryogenic heat exchanger was 
prepared under Contract NAS 9-12208. A controls investigation made using 
the computerized model was reported in NASA CR 115569. In that report a 
recommendation was made for a simple control , sensing Freon outlet temperature > 
and employing a constant speed reversible electric motor to actuate the H 2 
feed valve. In the absence of detailed knowledge of the thermal dynamics of 
the Shuttle Freon loop, characteristics for the control were selected on the 
basis of assumed values for: 

1. The maximum possible rate of change in Freon temperature 
at the heat exchanger inlet, and 

2. The effect of Freon temperature changes at the heat 
exchanger outlet on other elements in the system, which 
determine control specifications. 


14 



Hamilton U 

__ D4V|%«OM op UNiTCo ilHCOflfT COAPOA^TiChsi 

Standard ft® 


SVHSER 6180 


A - Thermal Stresses - Minimal 
B - Vibration Stresses - 

Requirement = 24 G's + 10^ cycles equivalent sine 

Mounted as shown, with full weight assumed to be supported by 
outer tube (1.0 O.D. , .028" wall, AAS062) gives factor of safety 
>1.5 for 10? cycles. 



^ Weight equal to 20 lbs . , including insulation. 

C - Shock - not a determining factor requirement on 20 G's 


L 


CONCLUSION - There are no foreseeable structural problems in the design 
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STRESS SUMMARY - CRYOGENIC H 2 HEAT EXCHANGER 
FIGURE 3 
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FUEL 

CELL 

HX's 

C3D 



CABIN WATER LOOP 
= 600 lb/hr 


FREON LOOP TEMPERATURES AND HEAT LOADS 
AT ASSUMED DESIGN POINT 

FIGURE 4 
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HYDRAULIC 
HX ^ 


(PLINE 5 - 41.6* 
■-18 lbs. 


CDLINE 6 - 110' 


CD LINE 3 - 87.3' 




FUEL 

CELL 

HX’s 


1“ 

“ “ 

=3 

r- 


“La 



n 

J 





□ 



i 


RAD BYPASS LINE - 20.8’ 


-29 lbs /unit 


GSE 

HX 


■13 lbs 


COLD ' 
WALLS 1 
334* LINE! 


CDLINE 2 - 37.4’ 


CDLINE 4 - 45.7’ 


OLINE 1 - 25’ 



OLINE C, CABIN WATER LOOP 
LINE TO COND. HX - 
70.8’ 

FREON LOOP LINES - 1" O.D. X .028” wall aluminum 

metal wt per ft = ,102 lb 

WATER LOOP LINES - 1” O.D. X .028" wall stainless st. 

metal wt per ft = .300 lb 


FREON LOOP LINE LENGTHS AND UNIT WEIGHTS 
FIGURE 5 
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TABLE I 

DYNAMIC PARAMETERS 


3 9 > 
3 3 

Q. — 

If 



NTU* 

FREON 

TRANSPORT 

TIME, 

Sec. 

OTHER 

FLUID 

TRANSPORT 

TIME, 

Sec. 

METAL 

HEAT**# 

CAPACITY, 

Btu/ a R 

A* 

FREON 

HEAT*## 

CAPACITY, 

Btu/°R 

OTHER 

FLUID 

HEAT*** 

CAPACITY, 

Btu/°R 

FREON 

hA, 

Btu 

HFTT 

OTHER 

FLUID 

hA, 

Btu 

FREON LOOP LINES, 
PER FOOT 

.0565 

.606 

- 

.0218 

.101 

- 

33.9 

- 

WATER LOOP LINES, 
PER FOOR 

.0687 

- 

1.82 

.0328 

- 

.304 

- 

41.2 

INTERFACE 
HEAT EXCHANGER 

10.60 

20.5 

61.9 

11.25 

3.43 

10.3 

6380. 

6380. 

FUEL CELL 
HEAT EXCHANGERS 
(3 IN PARALLEL) 

4.28 

23.4 

11.2 

4.41 

1.31 

1.87 

857. 

857. 

GSE HEAT EXCHANGER 

7.67 

4.1 

- 

2.78 

.68 

- 

4570 

- 

SUBLIMATOR 

4.43 

3.1 

- 

5.22 

.52 

- 

2660 

- 


* NTU = hA/WCp or UA/WCp in active heat exchangers 

** Half of metal plus stagnant fluid in active heat exchangers 

***Heat Capacity = Item Weight x Specific Heat = Q/AT 
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A maximum rate of change in Freon inlet temperature of 20°R per minute 
was assumed for normal operation, with possibly greater rates occurring 
during start-up if the radiator were suddenly by-passed. Hie control 
characteristic selected for. stability during normal operation was a full 
range of valve travel in 200 seconds and a + 1°R dead band. 

Under the present contract a mathematical model of the thermal 
dynamics of the Freon loop was prepared and programmed as an addition to the 
cyrogenic heat exchanger computer program. Using the expanded computer pro- 
gram the requirements for the cryogenic heat exchanger control were reexamined,, 
with particular attention to start -up conditions. 

The Freon loop system modeled was based on information obtained from 
the North American Shuttle ECS baseline (March 1972) and on studies conducted 
by Hamilton Standard, The active elements in the Freon loop during operation 
of the cryogenic heat exchanger with heat flows , fluid flow rates and temper- 
atures at the assumed design point are shown in figure 4. Line lengths and 
component weights are shown in figure 5. The model includes the Freon loop 
completed by the radiator bypass but not the radiator, hydraulic heat exchanger 
or associated lines. The line in the cabin water loop from the interface heat 
exchanger to the cabin condensing heat exchanger is also included in order to 
determine the effect of Freon loop temperature changes on the temperature of 
cooling water entering the condenser. 

The mathematical model and computer program are discussed in the next 
section and the controls reevaluation in the following one. 

Description of Math Model and Computer Program 

The components in the Freon loop are one of two types from the point 
of view of thermal dynamics : 

1. Those in which heat transfer occurs only between Freon 
and metal, i.e. transport lines and the inactive GSE 
heat exchanger and sublimator which have no second fluid, 
flowing. Freon and metal in these components come to the 
same temperature at steady state. 

2. Active heat exchangers (interface and fuel cell heat 
exchangers) . 

Component models for each type were selected for their simplicity and 
tolerance to input changes. Where required for accuracy, components were 
divided into two or more segments . 
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Nomenclature : 


A 

Cp 

m 

NTU 

Q 

T 

t 

At 

U 

W 

E 

@w 

T 

e 

u> 

j 


Heat transfer surface area 

Specific heat 

Mass 

I 

UA/W (Number of Thermal Units) 

Heat flow 

Temperature 

Time 

Finite interval in time 

Heat transfer coefficient (between two fluids 
in heat exchangers , between fluid and wall in 
lines) 

Flow rate 

Heat exchanger effectiveness 
Wall time const = (mCp)wall/UA 
Fluid transport time 
Base Napierian Logarithm (2.7183) 

Frequency (radians) 

pr 


Transport Lines and Inactive Heat Exchangers 

The fluid temperature at the outlet of a line with a uniform tube wall 
temperature is : 

T fluid out = T fluid e _NTU + T wall (l-e -1 ^™) 

The rate of change in average wall temperature , assuming no heat 
transfer from the outside of the tube, is: 

d T wa ll Q to wall (W Cp) fluid (T in - T out) 

3t (mCp)wall “ ~ ^ (nCp)wall 
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Employing the approximation that the assumed uniform wall temperature 
in the first equation changes as the average, these two equations may be 
combined to yield an equation for rate of change in fluid outlet temperature: 


fluid 
dT out 

— at — 


fluid , . . , 

dT in e _NTU , (WCp) fl uid (l-e _NTU ) (t 
3t (mCp)wall 


fluid_ T fluid \ 
in " 1 out / 


When using this equation, the fluid outlet temperature must be delayed 
by the time it takes a particle of fluid to pass through the line before 
it becomes the inlet temperature to the next element. 

The approximation involved in this equation for rate of change in out- 
let temperature does not introduce significant error when the NTU is on the 
order of 2 or less, as may be confirmed by calculating frequency response 
using this equation and comparing it to frequency response calculated using 
the exact transfer function for insulated lines. 


The above equation converted to a transfer function is : 


T fluid out 
T fluid in 


/ (WCp) fluid ^ 

= \ jqje'Nn 1 + (l-e~WTU) fmgpjwall J q-jwt 


iw 


(*- 





while the exact transfer function is : 


T fluid out __ 
Tfluid in 


jw 


NTU 

r + 9w~ 

Ct) 2 +1 

9w2 


+ NTU q? 

o>2 + 



Figure 6 shows a comparison of frequency response calculations using 
these two equations. Three lengths of Freon loop lines are shown. The 
NTU, which is proportional to line length, reaches a value of 2 at a line 
length about 36 feet. The approximate equation used in the computerized 
math model provides a reasonably accurate frequency response for a 30 foot 
line length, diverges somewhat from the exact solution at 50 feet, and 
becomes very inaccurate at a length of 100 feet. 

The logic in modeling the inactive heat exchangers in the Freon loop 
(sublimator and GSE heat exchanger) is precisely the same as for the lines, 
since the only significant heat transfer occurs between Freon and heat 
exchanger metal. Due to the large heat transfer surface area in heat 
exchangers, however, the NTU is typically large. Consequently, it is 
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FREON LOOP LINE TEMPERATURE TRANSFER FUNCTIONS 

FIGURE 6 
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necessary to break heat exchangers into segments having NIU's on the order 
of two. The sublimator was broken into two segments and the GSE heat exchang- 
er into four. Of the transport lines in the model, only the water loop 
line from the interface heat exchanger to the condenser required more than 
one segment. 

The differential equation for fluid leaving each line or inactive heat 
exchanger segment was put into the following form for solution by finite 
difference techniques using a computer : 


T fluid r-n = t 
out 1 


fluid 

out 


CD + 


fluid out 
T previous (n) 
segment 


_ t fluid ni 
out J 


x 


WCp) fl uid 
mCp) wall 


(l- e -NTUj x 


+ 


_ fluid out 
T previous 
segment 


. . _ fluid out 

(n-1) - T previous 
segment 


(n) 


x e'NTU 


where n = t segment 
At 


Values indexed n are values delayed by the fluid transport time in a 
segment, and are generated from currently calculated values (subscripted 
(1)) by setting: 


T (n) " T (n-1) 

for each of 2 through n values at each time step At. 

Active Heat Exchangers 

The steady state temperature of the cold fluid at the outlet of a 
counterflow heat exchanger, in terms of effectiveness, is: 

T cold „ E / T hot _ T cold\ + T cold 
out \ m m f m 

The steady state heat flow to the cold fluid is therefore : 

Qs s - WCp (t cold . T cold\ = WCp E (t hot - T cpld) 

\ s s s s / 
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The excess of the instantaneous heat flow to the cold fluid over the 
steady state heat flow at the instantaneous inlet temperatures determines 
the rate of change of heat exchanger wall material (one half of which is 
assigned to each fluid) : 

-CmCp) wall d T d ^ a11 = WCp (t cold _ x cold) . W Cp E ( T - T c g d ) 


d T wall _ fWCp) cold 
3t (mCp) wall 


,C1-E) 


T cold + P T hot _ T cold 
m m 1 out 


Writing, as for the transport lines: 


T cold _ x cold e -NTU + T wall (i-e _NTU ) 
out in \ f 


out in 

and combining this equation with the one above: 


d x cold d x cold , 

out _ m e _NTU + (i-e-MIu) 


~dt~ 


“3T 


(WCp) cold 
(mCpJwall 


(1-E) T c g d + E T - T c ° ld 


in 


out 


An equation for the rate of change of the hot fluid temperature leav- 
ing the heat exchanger is similarly developed. As with the equation for 
transport lines, the accuracy of this equation may be tested by comparing 
calculated frequency response to that calculated using an exact counterflow 
heat exchanger transfer function. An exact temperature transfer function 
was derived by Gyftopolous and Smets , and presented by Ball in "Approximate 
Models for Distributed- Parameter Heat Transfer Systems", ISA Transactions, 
Vol. 3, No. 1. A straight through transfer function for a heat exchanger 
with an NTU of 1, and fluid transport times and wall time constants equal 
to 10 seconds for both fluids has been replotted from Ball's paper in 
figures 7 and 8. 

The straight through transfer function resulting from the approximate 
equation is : 


t cold 
1 out 


T CQld 
m 
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NTU = 1.0, TRANSPORT TIME AND WALL TIME CONSTANTS 
FOR EACH FLUID » 10 SEC 



SOLID LINES - EXACT TRANSFER 
FUNCTION 

I 

SYMBOLS - COMPUTER MODEL 



0.01 0.1 

FREQUENCY - HERTZ 


HEAT EXCHANGER TEMPERATURE TRANSFER FUNCTIONS 
AMPLITUDE RATIO VERSUS FREQUENCY 


FIGURE 7 
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HEAT EXCHANGER TEMPERATURE TRANSFER FUNCTIONS 
PHASE LAG VERSUS FREQUENCY 

FIGURE 8 


26 



Hamilton 

Standard 


U 

ppvp&c*+ op ymTEo aircraft coetf*or*ATton 

R® 


SVHSER 6180 


Values calculated using this equation are plotted as syitibols in 
figures 7 and 8 and show good agreement with the exact frequency response.. 
The accuracy decreases as NTTJ becomes large, and high effectiveness heat 
exchangers were broken into several segments (the interface heat exchanger 
five, and the fuel cell heat exchangers two each). 

The finite difference form of the equation is: 



E T h <?t 
in 


,,, m tOlO, UUL 

(1) + (1-E) T previous 
segment 




out 


x At 


+ e'NTU x 


/_ cold, out cold out 

[ T previous (n-1) - T previous (n) 
\ segment segment 


and T^ = T^ n + ^ at each time step. 


Computer Program 


A listing of the computer program including the cryogenic heat 
exchanger and its controls, and the thermal dynamics of the Freon loop 
is given below. The input is defined by comment cards at the beginning 
of the listing. Values for the dynamic parameters in the Freon loop used 
as input to the computer program in this study are summarized in Table I. 
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TIME 

HOUR/MINUTES/SECONK 

09/21/09 


C DECK 505 TUBE IN TUBE COUNTEBFLOW HEAT EXCHANGER WITH SUPERCRITICAL 

C H2 ENTERING INNER TUBE » FREON IN ANNULUS. CONTROL 

C OF FREON OUTLET TEMPERATURE BY MODULATING H2 FLOW USING CONSTANT 

C SPEED REVERSIBLE MOTOR TO ACTUATE INLET VALVE 

c shutle freon loop thermal dynamics added 

c 

c input is loaded by standard load routine in the following locations 
c geometrical factors 

C 1 DI INNER DIAMETER OF INNER TUBE, INCHES 

C 2 TI THICKNESS OF INNER TURF WALL, INCHES 

C 3 DO INNER DIAMETER OF OUTER TUBE, INCHES 

C A TO THICKNESS OF CUTER TUBE WALL, INCHES 

C 5 XL HEAT EXCHANGER LENGTH, FEET 

C 6 CD COIL DIAMETER, INCHES 

C OPERATING CONDITIONS 

C 7 WD DESIGN VALUE OF H2 FLOW RATE, LB/HP 

C 8 WI INITIAL VALUE OF H2 FLOW RATE, LBS/HR 

C 9 WF FREON FLOW PATE, LB /HR 

f. 10 PD H2 PRESSURE AT HEAT EXCHANGFR INLET, PSIA 

C 11 ZNPTl NUMBER OF H? INLET TEMPERATURES VS TIME POINTS INPUT 

C 31-40 T1TIME TIME AT WHICH H2 INLET TEMP IS DESIGNATED, 1 THRU 10, SEC 

C 41-50 TITEMP H2 INLET TEMP AT TIME POINT 1 THRU 10, DEG R 

C MATERIAL PROPERTIES 

C 13 CPF FREON SPECIFIC HF AT , BTU/L8-0EG R 

C 14 CONDF FREON THERMAL CCNDUCT IV TTY, RTU/HR-FT-DEG R 

C 15 RHOF FREON DENSITY, LB/FT3 

C 16 VISE FREON VISCOSITY, LB/HR-FT 

C 17 CPW TUBE METAL SPECIFIC HEAT, BTU/LB-DEG R 

C 18 RHOW TUBE METAL DENSITY, LR/FT3 

C 19 XKW TUBE METAL THERMAL CONDUCTIVITY, BTU/HR-F T-DEG R , 

C n-P6 T3TEMP 18 INPUT TEMPERATURES FOR H2 ENTHALPY, DEG R 

C 91-106 T3FNTH H2 ENTHALPY AT DESIGN PRESSURE, BTU/LB 

C CONTROL VARIABLES 

C 20 GR MCTOR RATE, PERCENT VALVE TRAVEL PER SECOND 

C 21 SP FREON OUTLET TEMPERATURE CONTROL SET POINT, DEG R 

C 22 DB AND TEMP BAND AROUND SET POINT IN WHICH MOTOR IS INACTIVE, DEG R 

C 23 ACDL DELAY IN VALVE ACTUATOR FROM TIME OF TC SIGNAL, SEC 

C 24 TCTC THERMOCOUPLE TIME CONSTANT, SEC 

C COMPUTATIONAL CONSTANTS 

C 25 TM A X TIME AT WHICH TRANSIENT CALCULATION ENDS. SEC 

C 26 ZN NUMBER OF NODES IN HEAT F XC HANGER MODEL 

C FREON LOOP DYNAMICS 

C 111 FLHAF FREON LINE HA PEP FOOT , 8TU/ ( HR-FT-DEG R) 

C 112 FLMCWF FREON LINE HEAT CAPACITY PER FOOT, BTU/IET-DEG R1 

C 113 FLMCtF FREON HEAT CAPACITY IN FOOT OF LINE, BTU/1FT-DEG R1 

C 114 FL1L LENGTH OF FREON LINE 1, FT 

C 115 FL2L LENGTH OF FRFON LINE 2, FT 

C 116 FL4L LENGTH OF FREON LINE 4, FT 

C U7 FRP L LENGTH OF FREON RADIATOR BYPASS, FT 

C 118 STTRL LENGTH OF STAR TRACKER LINE, FT 

C 121 WLHAF WATER LINE HA PER FCOT, BTU / ( HR-FT-DEG PI 

C 122 WLMCWF WATER LINE HEAT CAPACITY PER FOOT, flTU/IFT— DEG R) 

C 123 WLMCLF WATER HEAT CAPACITY IN FOOT OF LINE, BTU/ 1 FT-DEG R1 

C 124 WLC L LENGTH OF WATER LINE FROM INTFR HX TO CONDENSOR, FT 

C 129 E I HX EFFECTIVENESS OF INTERFACE HX 

C 170 FIMCW WALL HEAT CAPACITY FOR FREON IN INTERFACE HX,BTU/OEG R 
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0 131 WIMCW WALL HEAT CAPACITY FOR WATER IN INTERFACE HX » BTU/DE G R 

C 132 TAUIF TRANSPORT TIME OF FREON IN INTER HX , SEC 

C 133 TAUIW TRANSPORT TIME OF WATER IN INTER HX, SEC 

C 134 EFC EFFECTIVENESS OF FUEL CELL HXS 

C 135 FCFMCW WALL HEAT CAP FOR FREON IN FUEL CELL HX$,BTU/DEG R 

C 136 FC4MCW WALL HEAT CAP FCR FC43 IN FUEL CELL HXS,BTU/DEG R 

C 137 TAUFCF TRANSPORT TIME CF FREON IN FUEL CELL HXS. SEC 

C 138 TAUFC4 TRANSPORT TIME OF FC A3 IN FUEL CFLL HXS, SEC 

C 139 GSENTU NTU OF GSE HX, 

C 140 GSEMCW HEAT CAPACITY OF GSE HX, BTU/OEG R 

C 141 TAUGSE TRANSPORT TIME OF FREON IN GSE HX, SEC 

C 142 SUB NTU NTU OF FREON I N SUBL IMATOR 

C 143 SU8MCW HEAT CAPACITY OF SUBLIMATOR, BTU/DEG R 

C 144 TAUSUB TRANSPORT TIME OF FREON IN SUBLIMATOR, SEC 

C 150 TIHXW INITIAL TEMPERATURE OF WATER INTO INTER HX, OEG R 

C 151 TFCWA INITIAL TEMPERATURE OF FC 43 INTO FUEL CELL HX A, OEG R 

C 152 TFCWBI INITIAL TEMPERATURE OF FC 43 INTO FUEL CELL HX B,DEG R 

C 153 TFCWCI INITIAL TEMPERATURE OF FC 43 INTO FUEL CELL HX C, DEG. R 

C 154 TBP! INITIAL TEMPERATURF OF RADIATOR BYPASS LINE, DEG R 

C 155 TGSEI INITIAL TEMPERATURE OF GSE HX, DEG R 

C 156 WW FLOW RATE IN WATER LOOP, LB/HR 

DIMENSION PL T f 50 2 , 1 3 1 

DIMENSION T 1 41 ) , TF(41), TWI41)., TSI411, ENT{ 41 ) 

DIMENSION 0(170), T1TEMPI10I, T1TIMEU0), T3ENTHI20), T3TEMPI 20 ) 
DIMENSION TL INI (91 ) ,TLIN2l 10 1 ) , TL I N4 ( 1 5 1) , TBP( 71 ) , TL INC! 501) 
DIMENSION TSTTRl 121 ) ,TIHXFO( 21 ), TIHXWI 181 ) ,TSU1 (21 ) , TSUB(2l) 
DIMENSION TIHXFK21 ) ,TI HXF 1 ( 2l),TIHXF2(2t) » TIHXF3 (21 ) ,TIHXF4(2t) 
DIMENSION TIHXWI (81 ) ,T IHXW2I 81 ), TIHXW31 81 ) , TIHXW4I 811, TIHXWO( 81) 
DIMENSION TFCFA1 (51) ,TFCFB1(51),TFCFC1 (51 ) , TFCFAO (51) , TFCFBO (51) 
DIMENSION TFCFCO( 51) ,TFCWA K 31 ), TFCWA1 ( 31) ,TFCWB1(31 ) ,TFCWCH3i) 
DIMENSION TGS1I25J, TGS2125), TGS3(25), TGSE(25) 

EQUIVALENCE 

X ( D ( 1 ) , DI ), t Dt 2 ) , T I ), ( D( 3 ) ,00 ), (D( 4), TO ), 

X . ( D l 51, XL ), ( D ( 6), CD 1, (0( 7*,WD ), ( D( 8),WI ), 

X (Dl 9 ) , W F ), (DUO), PD ), (D(ll>,ZNPTl ), (0(12>,ZNPT2 ), 

X 1D(13), CPF ), { D ( 14 ) , CONDF ), ( D( 15 ) ,R HOF ) , (0(16), VI SF ), 

X (D(17),CPW ), (Dl 18 ) »RHGW ), (D(19),XKW ), 10(20), GR ), 

X (0121), SP ), (0(22), DB AND ), (0(23), ACOL ), (0(24), TCTC ), 

X (0(25) ,TMAX ),(D(26),ZN ), (D (27), PLOT ), 

X l DC 31 > .TITIMEU)), 1 D (4 1) ,T 1TEMPU ) ) » 

X 1.0(71), T3TEMP11)), (0(91) » T3ENTHI 1 ) ) , 

X (D( 11 1 ) » FLHAF ) , (0(112), FLMCWF), (D(113), FLMCLF ), ( Dl 114),FL1L ) , 

X (0(115), FL2L ) , t Dl 1 lfc) ,FL4L ),(D(117),FBPL . ) , ( D( 1 18 ) ,STTRL ), 
X (0(121 ) , WLHAF ), (0(122), WLMCWF) ,(D( 123), WLMCLF) ,(D( 124) ,WLCL ), 

X <D( 129) ,EIHX ) , (D(130),FIMCW) , (D( 131) ,WIMCW), ( 0( 1 32 ) , T AU IF ), 

X (D( 133), TAUIW ) ,10( 134) ,EFC ) , ( DC 13 5) ,FCF MCW I , 1 0 (1 36 ) , FC4MC W ) , 
X ( 0(137), TAUFCF), (Dl 136 ) , T AUFC 4 ) , ( D( 139 > , GSENTU) , ( D (140 ) ,GSEMC W) , 
X (D( 141 ), TAUGSE) ,(D(142) ,SUBNTUI , ( D ( 143 ) , SUBMCW ) , ( D( 144 ) .TAUSUB ) , 
X (D(150), TIHXW ) ,(D( 151) , TFCWA I , ( D( 152) , TFCWBI ) , ( D( 1 53 ) ,Tf C WC I) , 
X ( D ( 1 54 ) , TBP I ) ,!P( 155), TGSEI ),(D(156),WW ) 

DO 12 1*1,170 
12 D( I ) = 0. 

15 CONTINUE 
‘ CALL START 
CALL LOAD(D) 


WRI TE l 6,999 ) DI, TI, DC, TO, XL, CD, WD, WI, WF, PD, CPF, CONDF, 
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X RHOF t VISFf CPW. RHOW» XKW, GRf SP, DBANDf ACDL , TCTC, FLHAF, 

X FIMCWF,FLMCLF,FL1L ,FL 2L * FLAL »FBPl , STTRL, WLHAF » WLMCWF , WLMCLF , WLCL , 
X E I HX , FIMCW, WIMCW, T All IF , TAUIW, EFC , FCFMCW, FC 4MC W, TAUFCF , 


X TAUFC4 

t GSENTU,GSEMCW, 

TAUGSE, SUBNTUt SUBMCW 

,T AUSUB, WW 



FORMAT ( 

*0 INPUT3/2X, 





XG12.5, * 

D I a,G12,5,» 

T I 3,G12.5,'DO 

3,G12. 5, ' 

TO 

a , 

XG12.5, * 

XL 3,G12 .5,* 

CD 3/2X, 




XG12.5, * 

WD 3,G12.5,* 

WI 3,G12.5,*WF 

3,G12.5,» 

PD 

3 t 

XG12.5, • 

CPF 3 f G12.5,» 

CONDF 3/ 2X , 




XG12.5,* 

RHOF 3,G12.5,* 

VISE 3,G12.5,*CPW 

3»G12,5,» 

RHOW 

8» 

XG12.5, < 

XKW 3 ( G12.5 f * 

GR 3/2X, 




XG12 .5, * 

SP 3,G12.5,» 

OBAND 3,G12.5, * ACOL 

3,G12.5,* 

TCTC 

3 1 

XG12.5, « 

FLHAF 3,G12.5,' 

FLMCWF3/2X, 




XG12.5, * 

FLMCLF3,G12.5, * 

FL 1L 3,G12.5,*FL2L 

3,G12.5,* 

FL4L 

a. 

XG 1 2 . 5 , ' 

FBPL 3,012.5,* 

STTRL 3/2X 




XG12.5, * 

WLHAF 3, G 12 . 5, * 

WL MC WF3 » G1 2 ,5, ' WLMCLF3 ,G1 2 . 5 , ' 

WLCL 

St 

XG12.5, * 

EIHX 3»G12.5, * 

FIMCW 3/2X 




XG12.5, » 

WIMCW 3,G12 .5, * 

TAUIF 3,G12.5,» T AU IW3, G12 . 5, * 

EFC 

a* 

XG12.5, * 

FCFMCW3, G12 • 5 , * 

FC4MCwa/2X 




XG12.5, * 

TAUFCF3, G12 .5, » 

TAUFC43,G12.5, • GSENTU3.G12. 5, 

' GSEMCW3 

XG12. 5, » 

TAUG SE3,G12.5,* 

SUBNTU3/ ?x 




XG12.5, ■ 

SU8MCW3,G12 .5, * 

TAUSUB3»G12*5» • WW 

a//) 



N=ZN 







IN 1 -N+ 1 
ITER=1 
I LP=0 
I PI =PLOT 
I IPL-IPL+1 
NPT1-ZNPT1 
VPOS=0. 

DVP0SA=0 .0 
TOFF =0. 0 
W=0. 

T I ME-0 • 0 

TFIIN1)=TGSEI 

ODI-D I ♦ 2*T I 

A - . 0054-5A+0 1 *D I 

AF =.005454*1 DO*DO— 00 l*OD I ) 

A$=.2618*DI*XL/N 

ASF=AS*0OI/DI 

AS$=.2618*DO*XL/N 

DH=D 1/12 

RHF = I 00- GDI )/ 12 

XMF=RHOF*AF*XL /N 

XMW=RH0W*.005454*(0CI*0D1-CI*0II*XL/N 
XMS=RHQW* ASS*TO/l 2 
DELT=XMF/WF*3600. 

TAUFLF^FLMCLF/ (WF*CPF1 *3600. 
FLNTUF=FLHAF/( WF*CPF) 

WLN TUF^ WL HAF/WW 
FNTUL1=EXP(-FLNTUF*FL U ) 

ENTUI H^EXPI -E IHX/I (l-EIHXl*5t ) 

ENTUL2= EXP ( -FLNTUF*FL 2L ) 

FNTUST =EXPt-FLNTUF*STTRL ) 

FNTUFC = I-EFC /( ( 1-EFC)*2) I 
FNTUBP= EXP ( - FLNTUF #FBPL 1 


30 



Hamilton U 

DWIAVUWTI4 t Co**0**t4** 

Standard p® 

IV G LEVEL 20.1 MAIN 


SVHSER 6180 

OATE = 72200 09/21/09 


ENTUL4* EXP I -FLNTUF*FL 4L 1 
ENTUGS = EX Pt -GSENTU/4 1 
ENTUSU=EXP(-SUBNTU/2) 

ENTULC = EXPI -WLNTUF*WLCL 1 
C LCX-FL MC WF*3600 • / (WF#CPF) 

CL I N 1=C LCC*F L lL/t 1-ENTUH) 

CLIN2=CLCC*FL2L/ t 1-ENTUL2) 

C $ TTR=C LCC* STTRL/1 1 -ENTUST 1 
CBP=CLCC*FBPL/( l-ENTUBP) 

CL I N4=CLCC*FL4L/ ( 1— ENTUL'4 ) 

CL I NC S WLMCWF* WLC L/ ( WW* 1 1 -ENTULC) >*3600. 

CIHXF=F IMCW/IWF*CPF*( l-FNTUI H 1 *5 1 *3600. 

C I HX W=W I MCH/ (WW*I1-ENTUIH>*5» *3600. 

CFCF=FCFMCW/( WF*CPF* ( 1-ENTUFC 1 *2 1 *3600 . 

CFC4=FC4MCW/ (WW*( 1-ENTUFC I *21*3600. 

CGSE*GSEMCW/ (4*WF*C PF* (1-ENTUGS ) 1*3600. 

CSUH = SUBMCW/(2*WF*CPF*U-ENTUSU1 1*3600. 
NIHXF=TAUiF/C5*DELTl 
NIHXW*TAUIW/15*DELT1 
NFCF=TAUFCF/(0ELT*2) 

NFC4=T AUFC4/ t OEL T*2 1 

NL I N1 = T AUFLF*FL1 L/OELT+N IHXF 

ML I N2=T AUFLF *FL2L/0EL T 

NSTTR-TAUFLF*$TTRL/DFLT+NFCF 

NBP=TAUFLF*FBPL/DELT 

NL IN4=T AUFLF*FL4L /OELT 

NLI NC=WLMCLF*WLCL/WW*3600./DELT 

NGSE=TAUGSE/t4*DELT) 

N$UB=TAU$UB/ (2*DELT 1 
TCEX=EXPIDELT/TCTCJ 
WTDR Y=( XMW + XMS 1 *N 
WTWET-WTDRY + XMF*N 
RW2=tOOI-OI> /( 48.*XKW) 

BFA-l.+3.5*( 00-0 0 1 1 /CO 
BFT=1, +3.5*01/00 
PRF = VI SF*CPF/CONDF 
REF=DHF*WF/(AF*VISF> 

I F (REF-2300 * 1 L . 1 . 2 

1 FFF=64. /REF 
HF=4.86*CONOF/OHF 

GO TO 3 

2 FFF=.184/(REF**.21 

HF= .023 *CONOF/OHF *REF **. 8*PRF**. 4*BFA 

3 DPF=FFF*XL/DHF*tWF/ AF) **?/ U .2 E+l 1 *RHOF 1 *BFA 
UFW= 1/1 1/HF + RW2) 

TCWF=XMW*CPW/ (UFW*ASF')*360C. 

TCSF*X MS*C PW/( UF W# A S S ) * 3 60 0 . 

TCFW=XHF*CPF/( UFW*ASF I *360 C. 

TCFS-X MF*CPF/ ( UFW*ASS 1 *3600. 

I F ( OE L T-TC WF1B5, 85,84 

84 WR ITE ( 6 » 1 000 1 TCWF, CELT 

1000 FORMAT { '0*** N IS TOO SMALL FOR STABLE SOLUTION ***3/5X ? 
X'TCWF#a,Gt3.5, 10X,'0ECT#3,G13.5) 

GO TO 15 

85 TOUT=550. 

TFOUT=500. 
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CALL UNLININPT1 »T1TIME,T1TEMP,TI PE, Til 1, I BUG 1 
CALL UN INTI 18 , T3TEMP ,T3ENTH » T( 1) , FNTI 1 ) , IBUG) 

DO 9 I = 1 » N 

T< I ■*■ 1 1 — T t I) + (TOUT-TU) ) /N 

TF ( I ) = TF I IN1 1 — €TF| INI l-TFOUT ) * IN- 1+1 ) /N 

TWt I 1=( TF( I 1+T 1 1 + 1) 1/2 

9 TS 1 1 ) = TF ( I 1 

10 TF I T -T F (1 ) 

TIT=TU Ml 

GO TO 21 

20 TF1-TFU1 
TFN1 = TF { I Nl 1 

CALL UNLIN(NPTI»T1TIME»T1TEMP»TIME»T(1) ♦ I BUG > 

CALL UN I NT U 8 , T3T EMP ,T3ENT H, T 111 , ENT t U , 1BUG1 

21 OP T0T=0 .0 

DO 100 1=1, N 
TFA B=( TFt I ) + TF (1*11 1/2 
IFU-O. 130,30,31 

30 DTW=(TFAB-TW IIII / TCWF 

GO TO 32 

31 TABMTtll+Tl 1+11 ) / 2 
TF I L= t T AB+TW (111/2 

IFt TFU-300. 141 ,41,40 

40 CPFIL=3 .55 
COND=.00196*TFIL**.637 

GO TO 42 

41 CPF IL= . 1805*TF IL**. 5223 
COND=.0002l*TFIL**l .028 

47 VI S=. 000263* TFIL**. 698 

PR = V IS *CP F IL /CONO 
PE=DH*W/(A*VI$1*TAB/TFIL 
H= .023*COND/nH*RE**,8*PP**.4*8FT 
UW= 1/(1 /H+RW2 1 
TC WH=X MW*C PW/ ( UW* AS | *3600. 

OENT=UW*AS/H* ( TWt I > -TAB 1 
ENT 1 I+l 1 = ENT ( 1 1 + DENT 

CALL UN1NTU8 , T3F NTH , T3T6MP, ENT (I + 1) ,T t T ♦ 1 1 , I BUG! 
TABMTt 1 1+TU + l) 1/2 
FF= . 1 84/ t RE **.2 1 
RH0=PD*144,/(55.1*TABI 

0P=FF*XL/(N*0H1*I W/Al**2/t 1. 2E+1 l*PHO) *BFT 
DPTOT=DPTOT+DP 

OTW = I TF AB-TWt I 1) /TC WF+ ( T A8-T W (I) l/TCWH 

32 TWt I ) = DTW*DELT/2+TWm 
DTS= ( TF AB— TS I I 1 ) /TCSF 
TS( I 1 =0 TS*DELT/2+TS (II 

OTF=(TW< I 1-TFABl/TCFW+tTSt Il-TFAB) /TCFS 
TF ( I )=DTF*DELT +TF( 1+1) 

1FIW-0. 133,33,34 

33 DTW= (TF AB-TWI I) l/TCWF 

GO TO 35 

34 DTW=(TFAB-TH( I 1 1 /TC WF+ ( TAB-TWtl 1 1 /TCWH 
3 5 TWt I ) = DTW*DELT/2+TW 1 1 1 

DTS=ITFAB-T$tl! l/TCSF 
100 TS( I ) = DTS*0£LT/2 + TS( II 

IF [TIME-DELT/2 1101,101, 110 
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101 IFCABStTFIT— TFI1) >-.01)102,102,103 

1C? I F < A 9 S t T I T-T I INI > )-.Cl > 1 04 * 104 , 1 03 

103 ITER= 1 

GO TO 10 

104 I TER= I TER + 1 

I F ( I T EP -N >10,10,105 

105 TCC =TF 1 1 1 
TF1=TFC1) 

00 777 LF= I , NL IN 1 

777 TL I N1 1 L F) =T F ( 1 > 

SEG IHX= El HX# ( TIHXW — TF1D1/5 
EIHX=EIHX/( EIHX+l l— E IHX ) *5 ) 

DO 782 LF = 1 , NI HXF 
TIHXFlt LFj^TFt 1)+SEGIHX*1 
TIHXF2 (LF> = TFm+SEGIHX*2 
TIHXF3I LF)=TFt 1 )+SEGlHX*3 
TIHXF4(LF)"TF( 1)+SEGIHX*4 

78? TIHXFO(LF)=TFt 1 ) +$EGIHX#5 

DO 787 LF =1 » NI HXW 
TIHXWI (LF)=TIHXW 
TIHXW1C LF1-T IHXW -SEGIHX*4 
TIHXW2I LF)=TIHXW -SEGIHX*3 
T IHXW3 ( LF )=T IHXW -SEGIHX*2 
TIHXW4(LF)=T IHXW -SEGIHX 

787 TIHXWOt LFl-TIHXW -SEGIHX*5 

DO 778 LF= 1 , NL IN2 

778 TLIN2C LF > =T I HXFO I N I HXF ) 

DO 786 LF= 1 , NSTTR 

786 T$TTR(LF>=TLIN2 INLIN2) 

SEGFCA=EFC*( TFCWA— TLIN2 CNLIN2) >/2 
SEGFC8 = EFC*I TFCWB I- TL I N2t NLI N2 ) )/2 
EFC=EFC/ ( EFCM 1-EFC >*2 ) 

DO 783 LF-i,NFCF 
TFCFA1(LF)=TLIN2(NLIN2>+SEGFCA 
TFCFB1 t LF ) =T LI N2 [ NL IN2 > +SEGFCB 
TFCFC1I LF)=TLIN2( NLIN2J+SEGFCB 
TFCFAO(LF)=TLIN2(NLIN2)+$EGFCA*2 
TFCFBOI LF > =TLI N2 ( NL I N2 ) +SEGFCB*? 

783 TFCFCOI LF)=TLIN2(NL IN2 ) + SEGFCB*? 
TFCFO=(TFCFAO(NFCF)+TfCFBO(NFCF)+TFCFCOINFCF) >/3 

DO 789 LF=i , NFC4 
TFCWAICLF)=TFCWA 
TFCWA1 ( L F >=T FCWA— SEGFCA 
TFCWB 1 ( LF ) =TFC WB I-SEGFCB 

789 TFCWC1(LF)=TFCWC!-SFGFCR 

DO 779 LF=1,NBP 

779 TBPILF>=TBPI 

DO 784 LF= l, NGSE 
TGS 1 1 LF )=TGS E 1 
TGS2ILF)=TCSEI 
TGS 3 CL F )=TGS El 

784 TG SE 1 LF ) = TGS El 

DO 780 LF = 1 , NLI N4 

78C TL I N4 ( L F) =TGSE I 

DO 785 LF “1 »N$UB 
TSU1(LF>=TGSEI 
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785 TSUBILF l = TGSEI 

DO 781 LF= 1 , NC INC 

7 81 TLINCt LE)=TIHXWGINIHXW1 

WR I TE 1 6 , 1 00 1 ) REF, HF , DPF, RE, H, DPTOT, WTDRY , WTWET, 

X TCWF, TCWH, TCSF, TCFW 
10C1 FORMAT ( *0INI TI AL COND I T I CN S3/2X, 

XG1 2 . 5, * REF 3, Gl 2 . 5 , * HF 3,G12.5,* DPF a,G12.5,' PE 3, 

XG1 2 . 5 , 1 H S,G12.5,* DPT0T3/2X, 

XG12.5,' WTD 3,G12.5,* WTW 3,G12.5»* TCWF3, G1 2 . 5 , * TCWH3, 

XG1 2 • 5, * TCSF 3,G12.5,» TCFW 3//) 

WRITEI6, 10031 

1 OC 3 FORMAT t • 0 H2 TEMPS, 1 THRU NS13I 
WR1TE(6,1002) IT (K»,K=1,IM) 

1 OC 2 FORMAT (2X,10G12.5) 

WR I TE ( 6 » 100A) 

ICO A FORMAT ( 'OWALL TEMPS, 1 THRU N3) 

WRI TE 16, 1002) lTWtK),K=l ,N) 

WRITEI6, 10051 

1CC5 FORMAT t *0 FREON TEMPS, 1 THRU N+13) 

WR! TE I 6,1002) < TF I K I ,K = 1 , 1 M ) 

110 TWF1 = TW( 1 )+UFW*RW2*lTF 1 1 l-TW(l) > 

PTLIN1= ITFt 1 l-TL INI (1)1 /CL INI 

TLINlt 11-TLlNlll )*OTLIM*DELT+lTFm-TFl)*ENTULl 
DO 888 LF= 2, NL TNI 

8 88 TLIN11 NUN1-LF+2>=TLINUNL INl-LF+1) 

DTIF1=( El HX* TIHXWI! 1 ) + U-E IH X) *TL I N1 (NL INI )-T IHX FI 1 1 ) ) / C IHXF*DEIT 
TIHXF1I 1 >=TIHXFltl)+DTIFl+(TLINl<NLINl — l)-TLINl(NLlNl) >*ENTUIH 
DTI F2= ( EIHX*TI HXW2 111 + 1 1-E IHX 1 *T IHXF1 1 N IHX F )-T I HXF2U) 1 /CIHXF*DELT 
TIHXF21 1 1 =T IHXF2 1 U + DTI F 2+ I T IHXF 1 1 NIHXF-l ) -T IHXF1 I N I HXF ) )*ENTU IH 
DTIF3= ( EIHX*TIHXW3{ 1)M 1-E IHX ) *T IHXF2I N IHXF )-T I HXF 31 1 ) ) /CIHXF*0ELT 
TIHXF31 l)=TI HXF3 1 1 ) +DTTF3+ (T I HXF2 ( NIHXF-l ) —T IHX F2 t N IHXF 1 )*ENTUIH 
OTIFA=( FIHX*TIHXWA< l)+( 1-E IH X ) *TI H XF3I N IHXF ) -T I HXF A ( 1 ) ) /CIHXF*DELT 
TIHXFAI 1 ) =T I HXF4 O ) *DT I F4+ t T I HXF 3 l N IHXF-1) -T 1HXF31 N IHXF ) )*ENTUIH 
DT I F5=( El HX* T I HXWI I 1)+ 1 1-E IHX ) *T I HX FA t N IHXF )-T IHXFOt 1 ) )/ClHXF*DELT 
TIHXFOl 1 )=T I HX FO ( 1 > +DT ) F 5+ ITIH XF 4! NIHXF-l ) - TIH XFAl NIHXF ) )*ENTUIH 
DTIW5=!EIHX*TIHXF4ll) + U- EIHX)*T IHXW II N IHXW )— T IHX W4( t M /CIHXW*DELT 
TIHXWAt 1)=TIHXW4( U+DTIW5+ t T I HXWI t NI HX W~U -TIHXWI 1 N I HXW ) ) *ENTU IH 
DT I WA= (EJHX*TIHXF3I 1 ) + ( 1-E IHX > *T IHXWAl N IHXW 1-TI HXW3 ( 1)1 /CIHXW*OFLT 
TIHXW3U) =TIHXW3I1 ) +DT 1 WA+ I T IHXWA ( N IHX W-l )-T IHXWAl N IHXW 1 >*ENTUIH 
DT I W3= ( E IHX* TI HXF2I 1 ) ♦ t 1-E IH X ) *UH XW3 ( NIHX H ) -T I HX W2 ID ) /Cl HXW* CELT 
TIHXW2I 1>=T IHXW2 1 1 ) +DT 1W3+(TIHXW3(NIHXW-1)— TIHXW3I NI H XW > ) *ENTUI H 
DTI W2=( EIHX*TIHXF1 (1) + U-E IHX ) *T I HXW? I NIHX W ) -T IHX W 1 ( 1 ) ) / C IHXW*DEL T 
TIHXWU U^TIHXWlll) +DTIW2+ 1TIHXW2I NI HX W-l > -T IHXW2 ( NI HXW )) *ENTU I H 
OTIWl={EIHX*TLINl<NLINl~NIHXF)+( l~E I HX ) *T I HXWI ( N IHX W )-T IHXWOtl) ) 

X /CIHXW*DELT 

TIHXWOI 1 ) = TIHXWOI 1J+DTIW 1+ ( T IH XW 1 1 NIHX W-l) - TIHX Wll N I HX W 1 )*E NTU I H 
TIHXWI (1)=T1HXW 
DO 903 LF= 2, NIHXF 

TIHXF1(NIHXF-LF+2)=TIHXF1INIHXF-LF+1) 

TIHXF2I NTHXF-LF*2)*T!HXF2t MHXF-LF + 1 ) 
TIHXF3(NIHXF-tF+2)=TIHXF3{NlHXF-lF+l) 

TIHXFAI NIHXF-LF*2)=TIHXFA(NIHXF-LF+ I) 

903 TIHXFOl NIHXF-LF+2 )=T1HXF01 MHXF-LF+1 ) 

DO 911 LF= 2* NI HXW 

TIHXWI (NIHX W-LF + 2 )=TIHXW1(M HXW-L F + 1 ) 

TIHXWAINIHXW— LF+2) =TIH XWA (MHXW-tF + 1) 
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TIHXH3I NIHXW-LF+2J*TIHXW3( MHXW-LF + 1 > 

T I HXW2 ( N I HXW-LF + 2 ) = TIHXW21 NIHXW-LF+1) 

TIHXWl ( NI HX W-LF+2 I =T I HXWlt MHXW-LF+1 ) 

911 TIHXWOt NlHXW-LF+2)=TIHXWG(MHXW-I.F+n 

DTL2DT* (TIHXFO(NIHXF)— TLIN2( 1)1 /CL IN2*DELT 

TLIN21 1 )=TU N2 tl) + DTL2 DT+ ( TI HXFO (N IHXF-1 )~T IHXFOI N IHXF ) ) #ENTUL 2 
DO 889 LF=2,NLIN2 

889 TLIN2I NLIN2-LF+2 ) = TL IN2INL IN2-LF + 1) 

DTSTTR=<TLIN2(NLIN2)~TSTTR11 M/CSTTR 

TSTTRI 1 >=T$TTR(1)+DTSTTR*DELT+(TLIN2(NLIN2-1)-TLIN2(NLIN2) )*ENTUST 
DO 913 LF=2, NSTTR 

913 TSTTRI NSTTR-LF+2 ) =T$TTP 1 NS TT R-LF + 1 ) 

DTFA1=(EFC*TFCHAII l) + t 1-EFC ) *T STTR ( NSTTR) -TFCF A 1 ( 1 M/CFCF*DELT 
TFCFA1 I 1 > -TFCFA1 ( 1 ) + DT FA 1 + f T STTR (NSTTR- 1 )~TSTTR ( NSTTR ) ) 4FNTUFC 
DTFB 1® ( EFC*TFC WB11 1 ) + ( 1-EFC)*TSTTP (NSTTR >-TFCFBl( 1 ) )/CFCF*DELT 
TFCFB1 ( 1 ) =TF CFfii 1 1) +DT FB1 + (TSTTR(N STTR-1 )- TSTTR ( NS TTR ) ) 4ENTUFC 
DTFC1=( FFC*TFCWC1 ( l )+( 1-EFC) *TSTTR (NSTTR )-TFCFCl( 1 ) >/CFCF*DELT 
TFCFCK 1)®TFCFC1( 1 ) +OTFC 1+ I T STTR ( N STTR- 1 ) -TSTTR ( NS TTR )) *ENTUFC 
DTFA0=( EFC*TFCWAI(1)+1 1-EFC > *T FC FA 1 ( NFCF )-TFCFAO( 1) ) /CFCF*DELT 
TFCF AOtl) =TFCFAOU) +OTFAG+ (TFCFA1 ( NFCF-1 I-TFCFAH NFCF ) )*ENTUFC 
DT F BO® ( EFC*TFC WB I + 1 l-EFC ) *TFCFB 1 ( NFCF I -TFCF BO f 1 ) ) /CFCF*DELT 

TFCFBO(l)®TFCFBO(l) +DTFBO+ (TFCFB 1 ( NFCF- l)-TFCFBl(NFCF)) *ENTUFC 
DTFCO=( EFC*TFC WC I + (1-EF C) *TFC FC 1( NFCF) -TFCFCOI l ) ) /CFCF*DELT 

TFCFCOI 1 l=TFCFCO I 1 > +DTFCO+ 1TFCFC t ( NFCF-1 )- TFCFC 1 1 NFCF) ) +FNTUFC 
DTWAl* ( EFC*TFCFA1 ( 1 )'+'( 1— EFC ) *T FCW A I (NFC A )— TFCWA H l ).) /CFC4*DELT 
TFCWA 1( 1) =TFCWA1 ( 1) +DTWA1+ IT FC WA It NFC4-1 ) — TFCWA ! ! NFC 4 )) * ENTU FC 
PTWB1= ( FFC*TFCFBH l)+( 1-EFC ) *TFCWB I -TFCWBl(l) )/CFC4*DELT 

TFCWBK 1)*TFCWB1 (1M-DTWB1 

DTWC 1= ( EFC*TFCFC 1 ( l ) +1 1-EF C) *TFC WC I -TFCWCl (1 ) )/CFC4*0ELT 

TFCWC1 ( 1)=TFCWC1 ( U+DTWCl 
TFCWAI ( 1) =TFCWA 
TFCFOO* TFCFO 

TFCFO* ( TFCFAO (NFCF ) +TFCFRO (NFCF) +TFCFCOI NFCF) ) /3 
DO 893 LF=2 • NFCF 

TFCFA1 ( NFCF-LF+2 )= TFCF A l (NFCF-LF+ 1 ) 

TFCFB1 ( NFCF-LF+2 )=TFCF 0 11 NFC F-LF + 1 ) 

TFCFCll NFCF-LF+2) =TFCFC1 (NFCF-LF+1) 

TFCFAO ( NFCF-LF + 2 ).=TFCF AO (NFCF-LF+1 ) 

TFCFBOI NFCF-LF+2 )-TFCFEOt NFCF-LF+1 ) 

893 TFCFCOI NFCF-LF+2) ®TFCF CO (NFCF-LF+1 ) 

DO 910 LF®2,NFC4 

TFCWAI ( NFC4-LF+2 1=TFCWAI(NFC4-LF+1 ) 

TFCWA 1 1 NFC4-LF + 2 ) = TFCWAI ( NFC4-LF+ 1 ) 

TFCWBl (NFC4-LF+2 )=TFCWB1 (NFC4-LF+1 ) 

. 9 1C TFCWC1(NFC4-LF+2)*TFCWC) INFC4-LF+1) 

DTBP=< TFCFO -TBP1 1) )/CBP 

TBP 1 1 ) =TBP( 1 )+DT BP*DELT + l TFCFO-TFCFOO) *ENTUBP 
00 890 LF=2»NBP 

89C TBP 1 NBP-LF + 2 )*TBP ( NBP-LF + 1 ) 

DTGS1 =( TBP( NBP)-TGS1 (1 M/CGSE 

TGS 1 ( l ) =TGS 1 ( 1)+0TGS1*DELT+(TBP( NB P-1 ) -TBP ( NBP ) )*ENTUGS 
0TGS2=( TGSt (NGSE )-TGS2ll) ) /CGSE 

TGS211) =TGS2in+0TG$?*DELT + ( TGS1 (NGSF-1 )-TGS H NGSE V ) *ENTUGS 
DTGS3* ( TGS2 ( NG $E ) - TGS3 ( 1) ) /CGSE 

TG$3U)*TGS3tl ) +DT GS3*DELT ♦( TG $2 (NGSE- 1 )-TGS ?( NGSF ) )*ENTUGS 
DTG SE* ( TGS3 ( NGSE )-TGSE ( 1 ) )/CGSE 
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TGSE! 1I=TGSE(1 M-DTGSE*DELTM TGS3 (NGSE- 1 l-TGS 3! NGSE ) > *ENTUGS 

DO 898 LF=2» NGSE 

TGSl (NGSE-LF+2 ) = TGS 1 ( NGSE-LF +1 ) 

TG S2 1 NG SE-LF+2 ) = TGS2 ( NGSE- LF +1 ) 

TGS 31 NGSE-LF + 2 ) = TGS 3{ NGSE-LF + 1 ) 

898 TGSE! NGSE-LF +2 ) =TG$ E! NGS E-LF + 1 ) 

DTLIN4=t TGSE! NGSE) -TLIN411 )) /CLIN4 

TLIN4( 1 )=TL I N4 1 1 ) +DTI I N4*DEL T+ (TGSE (NGSE- 1 )-TG SE t NGSE >) *ENTUL4 
On 891 LF =2 » NL T N4 

891 TLIN4!NLIN4-LF+2)“TLIN4!NL IN4-LF+1 > 

DTSUi-ITL IN4!NLIN4)-TSUIIll)/CSU8 

TSUI (1) =TSU 1 m+DT$Ul*DELT*!TL IN4INIIN4-1 )—TL I N41 NL I N4 ) )*ENTUSU 
DTSUB*( TSUKNSUBI-TSUBU) 1/CSUB 

TSUB! I ) =T SUB ! 1 >*DTSUB*DELT+< TSUI! NSUB- 1J-TSU1! NSUB) KENTUSU 
00 895 LF = 2 » NSUB 
TSUUNSUB-LF+2 ) = TSU 1( N SUB-LF+ 1 1 
895 TSUBI NSUR-LF+2 > =T$U B ( NSUB-LF ♦ 1) 

TF( INI ) -TSUB t NSUB I 

OTLCDT= (TIHXWO IN IHX W )-TL INC! 1 ) ) /CL INC*DELT 

TLINCt 1»=TLI NC 1 1 ) +DTLCDT+ I T I HX WO (N IHXW-1 )— T IHXWOI N I HXW ) )*ENTULC 
DO 892 LF=2t NL1NC 

892 TLINCI NL1NC-LF+2)“TLINC(NL 1NC-LF+1) 

IF II PL-IIPLU13 t 112,112 

113 1 1 PL= 1 

WRI TE 16 » 1007 I 

1007 FORMAT! '0 SECGNDS F2 FLOW TUI TIN+1) 

X TF( 1) TF ! N£1 ) TW (1) TWF1 3) 

WRITE(6»1 006 ) T IME* W t Till* T(INl), TF (1) , TF ! INI) * TWl 1) , THF 1 
WR I TE ! 6 » 1 008 ) 

1008 FORMAT I *0 INT HX FREON SEG 1 TO 4 

X INT HX WATER SEGS 4 TO 1 a) 

WRITE! 6 t 1006) T1HXF1 I NIHXF) ,T I HXF2 1 NIHXF ) , T I HXF3 (N IHXF ) , 

X T I HXF4 !N IHX F ) » T IHX W4! N IHXW) , TIHXW3! N I HXW ) , TIHX W2 ! N I HXW) , 

X TIHXW1 ( N IHX W) 

WRITE (6, 1006 I TSUI (NSUB) , TL IN4 l NL I N4 1 tTGSE (NGSE )♦ TGSl (NGSE) t 
X TBP INBP) ,TFCFAO(NFCF) .TFCFA1INFCF > , TSTTR (NSTTR ) 

WRITE! 6»1006)TFCWAI f NF C4 ) , TFCW At IN FC4 ) , OT FA 1 » DT WA 1 , T FCF B l! NFCF ) * 

X TFCWBKNFC4) » TFCFCGI NFCF ) 

ILP = ILP+1 

IF ( I LP .GE. 501) GO TO 112 
PLT! IL P » 1 ) = T IME 
PLT ! I L P » 2 ) =T L I N4 ! NL IN4 ) 

PLT ( I LP »3 ) = TF ! I N1 ) 

PLT! ILP ? 4 )- W 

PLTULP*5)*TFm 

PLT! ILP»6)=T LINK NL INI) 

PLT! iLP»7)-TIHXF0(NIHXF) 

PLT(ILPt8)=TSTTR! NSTTR) 

PLTULP,9)=TFCF0 
PLT! ILP » 10 )-T BP ! NBP ) 

PLTULP,11I=TGSE!NGSF) 

PLT! lLPt 12)=TIHXW0I NIHXW) 

PLT t ILP,13)=TL INC! NL INC ) 

112 I IPL*I IPL*1 
10C6 FORMAT! 2X.8G15. 5) 

TT- t TF1 +TF ( 1 ) ) /2 
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TCC = TT+ (TCC-TT J/TCEX 
TI ME=T I ME+DELT 

IFtTIME-TMAX >120, 120,200 
120 DVPOSO=DVPOS 

417 I Ft I TLIN4I NL I N4— 7 1 ) -TL I N4( NL IN4-72 ))/DELT+.4)418,418»419 

419 IF t t TL I N4t NL IN4— 7 1) - TL IN 4 ( NL IN4-72 ) I / DEL T-. 4)4171,421, 421 

4171 I F ( TCC-S P-DBAND/2 )4191,418,418 

418 DVPOS=GR 

GO TO 422 

4191 I F ( TCC-SP + DBAND/2 ) 42 l , 420, 420 
42C DVPOS-O.O 

GO TO 422 

421 DVPOS=-GR 

422 l F( D VPOSO-DVPOS ) 424,423,424 

423 I F ( T I ME-TOFF-ACDL 1426,426,425 

424 TOFF = T I ME 

GO TO 426 

425 DVPOSA=DVPOS 

426 VPOS=VPO$*DVPOSA*DELT 

IFIVPOS 1427,427,428 

427 VPOS=O.G 

GO TO 430 

428 IFtVPOS-1. 2)430,430,429 

429 VP0S=1 .2 

43C W=WO/. 966* SI Nt 1 . 57* VPO S/1 . 2) 

GO TO 20 
2 GO CONTINUE 

I F ( IL P .GE . 501) I L P=50C 
CALL CALULMIPLT, ILP,12,502) 

GO TO 15 
END 
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Controls Reevaluation 

The H 2 feed control selected for the cryogenic heat exchanger in a 
study conducted under Contract NAS 9-12208 had the characteristics indicated 
in figure 9: a constant valve stroke travel rate equal to 200 seconds for 

full valve travel to the design flow rate, and a two degree dead band around 
the control set point. The H 2 source considered for the cryogenic heat 
exchanger was a supercritical tank for fuel cells. H 2 is tapped off an 
external pressurization system for this tank as shown in figure 10. The Freon 
temperature leaving the cryogenic heat exchanger is sensed, and the control 
is activated when the sensed temperature is one degree more or less than the 
500°R set point. 

Conservative non-ideal control characteristics assumed for the controls 
study were: 

Thermocouple lag ■ 1.0 second 

Actuator delay =0.2 second 

Non - linear flow versus valve stroke: 

W 1 / 1.57 _ . \ 

W design = 79 ^ (Sin 1 7T x stroke j 

The non-linear Hj flow versus valve stroke may arise due to valve 
design or from variations in pressure in the supply tank which changes 
flow versus stroke relationships by changing the density of H? passing 
through the valve. With a reasonably careful valve design, the non-linearities 
due to pressure changes should be within those assumed because the pressure 
in the relatively large supply tank cannot change at a high rate. 

The type of control described was selected for its simplicity, since 
it requires only a constant speed reversible motor for its actuator; and 
for its low power consumption, since it requires no actuator power at steady 
state conditions. Using this type of control, the valve travel rate must be 
restricted to full travel in 200 seconds for stable operation given the tube - 
in - tube heat exchanger dynamic characteristics and the assumed non-ideal 
control characteristics. However, this valve travel rate was found to be 
adequate to control Freon outlet temperature with rates of change in Freon 
inlet temperature up to 20°R per minute. A 20°R per minute ramp was thought to 
be more severe than could occur in the Freon loop except during startup. 

Further study was recommended for the startup situation. 

Under the present contract, the computerized model of cryogenic heat 
exchanger, controls, and Freon loop thermal dynamics described in the pre- 
vious section was used as a tool to reevaluate the assumptions made in the 
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control selection, and to provide a definition of requirements during 
startup. Specific goals for this investigation were to: 

• Confirm that during normal operation the rate of Freon temperature 
change at the cryogenic heat exchanger inlet does not exceed 20°R 
per minute for any appreciable time, and consequently that the sel- 
ected control provides adequate response during normal operation. 

• Determine whether Freon exit temperature perturbations caused by 
control action can return around the Freon loop in such a way as 
to be detrimental to control stability. 

• Determine the effect of Freon exit temperature perturbations on 
the temperature of cooling water entering the cabin condensing 
heat exchanger and Freon entering the fuel cell heat exchangers. 

This determination allows a more realistic specification of control 
limits, especially during startup. 

• Define the startup transient and modify the control as required 

to provide acceptable temperature control limits during the startup. 

TVo input disturbances were applied to evaluate the action of the 
cryogenic heat exchanger control under non-startup conditions: 1) starting 

with design heat loads, one fuel cell was "shutdown" by lowering the 
temperature of FC 43 W entering one fuel cell heat exchanger in a two min- 
ute period; and 2) again starting with design heat loads, one fuel cell 
and the cabin water loop were simultaneously "shutdown" in a two minute 
period. The second disturbance is more severe than can actually occur 
because the thermal mass of the water loop will prevent cooling of the 
water at this high rate. 

The startup transient was investigated by setting design heat loads 
on the interface and fuel cell heat exchangers, with the Freon in the loop 
from the beginning of the radiator -bypass through the inlet to the inter- 
face heat exchanger at S00°R. This is the situation at the instant the 
radiator bypass valve is opened. 

The results of these investigations are presented in computerized plots 
in figures 11 through 14. The location of each variable plotted is located 
by number on the schematic in figure 5. Some general remarks about certain 
features of the thermal dynamics of the system may help in interpreting 
these plots: 

• Transport time of the liquid around the loop is on the order of 
150 seconds. Heat added to a particle of Freon cannot result in 
a temperature rise in a downstream component until after the 
period of time it takes the fluid particle to travel to that 
component. 

nr FC 43 is a 3M Company coolant. 
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TRANSIENT RESPONSE - SHUTDOWN ONE FUEL CELL 
FIGURE 11-A 
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TRANSIENT RESPONSE - SHUTDOWN ONE FUEL CELL 

FIGURE 11 -B 
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TRANSIENT RESPONSE - SHUTDOWN ONE FUEL CELL 


FIGURE 11 -C 
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TRANSIENT RESPONSE - SHUTDOWN FUEL CELL 8 H 2 0 LOOP 


FIGURE 12 -A 



















TRANSIENT RESPONSE - SHUTDOWN FUEL CELL § H 2 0 LOOP 

FIGURE 12 -C 
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TRANSIENT RESPONSE - STARTUP WITHOUT LEAD 
FIGURE 13-A 
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TRANSIENT RESPONSE - STARTUP WITHOUT LEAD 
FIGURE 13-B 
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TRANSWITH RESPONSE - STARTUP - MODERATE UEAD 
FIGURE 14-B 
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• The fuel cell heat exchangers and especially the interface heat 
exchanger are highly effective. The steady state change in the 
temperature leaving a heat exchanger is one minus the effective- 
ness times the change in incoming temperature of that fluid 
(assuming no temperature change in the other fluid) . This 
value is .085 for the interface heat exchanger. 

• The sum of thermal mass in the Freon loop is large, and the sum 
of metal time constants around the loop is of the same order of 
magnitude as the control valve travel time. 

Each of the plots is discussed below. 

Shutdown of One Fuel Cell (Figure 11) 

Starting with the design heat load on the cryogenic heat exchanger, 
the shutdown of one fuel cell was simulated by linearly reducing the 
temperature of FC 43 coolant entering one fuel cell heat exchanger from 
629. 2 “R to 553. 2 “Rover a period of two minutes beginning at time zero. 

The resulting temperature transient and control action are shown in figure 


The maximum rate of change of Freon temperature at the cryogenic heat 
exchanger inlet is about 8“R per minute, well under the 20“ R per minute the 
control was designed to handle. Consequently the Freon outlet temperature 
does not go outside the 499“R to 50i°R control dead band by more than 0.3°R, 
well within tolerances of ± 5“R. 

A comparison of Freon tenperatures entering and leaving the interface 
heat exchanger (plots 8 and 9) shows that none of the "high frequency" 
temperature fluctuations caused by control action are able to pass through 
the large, high effectiveness interface heat exchanger. Therefore, there 
is no possibility of unstable control action resulting from temperature 
feedback around the Freon loop. 

Simultaneous Shutdown of One Fuel Cell and Cabin Water Loop (Figure 12) 

Figure 12 shows the results of reducing the inlet FC 43 temperature of 
one fuel cell from 629. 2 “R to 553. 2 “Rand reducing the water temperature into 
the interface heat exchanger from558.2°R to 500 “Rover two minutes starting 
at time zero. Although the transient lasts much longer than when a single 
fuel cell is shutdown, the rate of change in Freon temperature at the 
cryogenic heat exchanger inlet is scarcely greater. This is the result 
of the transport delay between the time a particle of Freon enters the 
interface heat exchanger and the time the same particle enters the fuel cell 
heat exchanger. 
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Only the initial portion of the plotted transient (betwen 300 and 400 
seconds) is a realistic representation of the response to the input described 
above. This is so because the fuel cell cooling loop has not been included 
in the computer model and the FC 43 temperature at the inlet to the active 
fuel cell heat exchangers was constant. Without adjustment in FC 43 loop 
temperatures to changes in Freon temperature, the correct steady state is not 
approached at the end of the transient. 

The input disturbance in .this case is more severe than can be imagined 
to occur in reality, because the large thermal mass in the cabin water loop 
would prevent water temperature reduction at the rate input. Consequently 
the control described at the beginning of this section is more than adequate 
in non- startup situations. 

Startup (Figures 13 and 14) 

The most severe Freon. temperature changes at the cryogenic heat exchanger 
inlet will occur when the radiator is effectively removed from the Freon loop 
prior to descent, and the. full heat rejection load is. imposed on the hereto- 
fore inactive cryogenic heat exchanger. 

The cryogenic heat exchanger H 2 control valve will initially be closed 
and will require 200 seconds to get to a full-open position. 

The radiator may be removed from the system by immediately bypassing 
Freon around the radiator, or by first stowing the radiator and bypassing it 
at some later time. If the radiator is stowed before bypassing, the heat load 
will be inposed on the cryogenic heat exchanger more gradually because the large 
radiator metal mass will provide cooling for a period of time after heat rejection 
by radiation to space has stopped. At the present time, there appears to be no 
disadvantage to first stowing the radiator. However, since bringing the 
cryogenic heat exchanger on-line by suddenly bypassing the radiator is the most 
severe startup mode, it is the one considered in this study. 

The startup transient that results with use of the control previously 
described is shown in figure 13. The Freon tenperature at the cryogenic heat 
exchanger exit rises to 545°R, or 45® above the control temperature. It is 
outside the control dead band for a. period of about 200 seconds. The effect 
of this startup temperature excursion an the operation of the fuel cells and the 
the cabin condensing heat exchanger may be; estimated from figure 13. Figure 
13 shows the Freon temperature at the inlet to the fuel cell heat exchangers 
rises less than a degree, and therefore, startup will have negligible effect 
on the fuel cells. The small tenperature rise at the fuel cell heat exchanger's 
inlet occurs because the interface heat exchanger acts as a buffer. 
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The cooling water entering the cabin condensing heat exchanger, shown 
in figure 15, rises 30°R above its design value to 53S°R. This will cause a 
temporary rise in cabin humidity because some of the water which has been 
condensed in the heat exchanger will be re-evaporated. The duration of 
high cooling water temperature is so short, however, that it is unlilcely 
that serious discomfort will result. 

Although the startup transient shown in figure 13 may be acceptable, 
the magnitude of temperature excursions during startup can be reduced by 
the addition of simple anticipation or lead to the controller. By obtain- 
ing the time derivative of Freon temperature at some point upstream of the 
cryogenic heat exchanger, appropriate control action can be taken before 
a temperature change is sensed at the cryogenic heat exchanger Freon loop. 

As an example, two theraocouples could be placed a measured distance apart 
in line 4 (see figure 5) . The difference between the two thermocouple 
signals would be proportional to rate of change in Freon temperature which 
will appear at the cryogenic heat exchanger inlet at some later time. When 
the rate of change exceeds a certain value, the controller causes the actuator 
to begin opening the H 2 control valve. This anticipation is required only 
during startup so the rate of change that activates the control would be 
set at 20°R per minute. 

Figure 14 shows the startup transient using such anticipation in the 
controller. The temperature rise in cooling water at the condenser inlet 
is reduced to 19 degrees, with a maximum temperature of 524°R occurring. 

The locations of the anticipation thermocouples in this case were about 5 
and 5-1/2 feet upstream of the sublimator. These locations appear reasonably 
close to optimum since an attempt to increase the amount of lead by moving 
the thermocouples 29 feet upstream produced the control action plotted in 
figure 15. The valve started to open earlier (as indicated by the H 2 flow 
rate) , but the temperature derivative went below its threshold value before 
the temperature at the cryogenic heat exchanger exit rose to the lower limit 
of the control band. As a result, the valve started to close again, the 
exit temperature rose sharply, and the resulting temperature excursion is 
of the same magnitude as in the previous case. 

Conclusions of Study 

• A constant speed control valve actuator with a rate of full travel 

in 200 seconds provides adequate temperature control in all non-startup 
situations. 

• During the most severe startup mode such a control allows a 30 <> R 
temperature rise in cooling water entering the cabin condensing 
heat exchanger. The temperature is above its design value for a 

period of 200 seconds. This short time excursion seems acceptable, 
but its magnitude may be reduced to 19^ or less by adding anticipa- 
tion to the controller. 
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TRANSIENT RESPONSE - STARTUP - MAXIMUM LEAD 
FIGURE 15 -B 
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• No temperature fluctuations caused by control action can feed back 
through the Freon loop. 


Heat Exchanger Breadboard Design 

Using information and stress analysis conducted under the design study, 
the heat exchanger was depicted in drawing SVSK 86021. This design was evolved 
after discussions with Jones Metal Products Company, who fabricated the final 
heat exchanger. The design was reviewed with the NASA on August 23, 1972 and 
approval to manufacture was obtained. 

Heat Exchanger Procurement 

Requests for bid on the tube - in - tube cryogenic heat exchanger depicted 
on SVSK 86021 were sent to three outside vendors and to the Hamilton Standard 
manufacturing facility. Only one affirmative reply was received, from the 
Jones Mstal Products Company of East Windsor, Connecticut. The remainder 
either did not have the specific tube rolling capability or did not have the 
current shop capacity to handle this job. Therefore, a purchase order was 
placed with Jones Metal Products Company . 

Unfortunately twelve feet was the maximum tubing length available without 
purchasing a mill run. The first attempt at manufacture was to butt weld two 
lengths of tubing. The initial attenpt to manufacture the heat exchanger 
disclosed two problems: 

1. Fracture occurred during coiling at the butt weld used to extend the 
tubing length. 

2. Severe rippling of the outer tube occurred at the inside of the coil. 

In order to eliminate the weld fracture problem the heat exchanger length 
was reduced as shown in SVSK 86021, Revision A. As can be seen in figure 2 , 
this results in a configuration with a system weight approximately one pound 1 2 * * 5 
more than optimum but allows the basic objective of the program, which is cor- 
relation of theimal and dynamic analysis with test. Should prototype implementa 

tion be required, time and funds should be alloted to obtain a mill run of the 
proper length. 

To correct the rippling problem the thickness of the outer tube was 
increased from .028 inches to .035 inches. This resulted in a marginally 
acceptable condition where minor dents in the tubing served as the focal 
point for rippling. Fortunately with careful rolling and working of the 
tubing rippling was held to an acceptable limit. Future designs utilizing 
this tubing should have the coil diameter increased from 10.5 inches to about 
12.0 inches. 
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The heat exchanger was received at Hamilton Standard and subjected to 
300 psig leakage/pressure tests and was found acceptable. The completed 
heat exchanger is depicted in figure 16 and the x-ray shown in figure 17 
shows the separation between inner and outer tubes. It was then cleaned 
to HS 3150 CE-3 in the hydrogen loop and HS 2241 G in the Freon loop, 

Set-Up Test Facility 

Upon completion of the system dynamic analysis, meetings were held with 
design, analysis and facilities engineering to define the detailed test require- 
ments and requisite facilities. The test requirements are contained in 
.Appendix A. It should be noted that any valves required to perform transients 
were considered as part of the test equipment and no simulation of the closed 
loop temperature control valve, controller, or sensor would occur within this 
test program. The transients defined in Appendix A are acceptable to evaluate 
dynamic response and obtain a thermal transfer function. 

It became apparent from review of the detailed test requirements that test 
costs would greatly exceed original estimates. Since the major element of these 
costs was the test set-up, it was decided to investigate test houses who already 
had these facilities. Contact was made with Linde, Inc., Wyle Laboratories, 

Pratt § Whitney Aircraft of Florida, Ogden Laboratories, and Beech Aircraft of 
Boulder, Colorado. Firm bids were received from Beech Aircraft and Wyle Laboratories. 

Wyle Laboratories provided the lowest bid as well as taking no exception 
to the specification. Beech Aircraft wanted to use Freon 11 in place of 
Freon 21 and wanted subcritical rather than supercritical hydrogen. Consequently, 
a purchase order for the test program was placed with Wyle Laboratories. 

Wyle prepared a test procedure which is included in the Master Test Plan 
contained in Appendix B. The test plan was reviewed by NASA on October 11, 1972 
and a final plan including NASA comments was submitted on October 24, 1972. 


Test Set-Up arid Instrumentation 

The test set-up is shown schematically in figure 18. The test heat 
exchanger (specimen) and hydrogen handling equipment were located out-of- 
doors as in figures 19 through 22. Hydrogen was obtained in subcritical 
liquid form and contained in the 2000 gallon vacuum insulated tank shown in * 
figure 20. Boil off from the tank entered a 5 stage compressor, was com- 
pressed to approximately 2300 psig and contained in cascade tanks. Pressure 
was regulated to 500 psig in an ullage tank and regulated to a value above 
the critical pressure prior to entering a cooling coil located in the 2000 
gallon tank. Cooled supercritical hydrogen exited the coil and was directed 
to the specimen inlet or a dunrny specimen while not operating the specimen. 
Hydrogen inlet and outlet pressures were measured on 0-300 psi gages. 
Hydrogen inlet and outlet temperatures and differential pressure across the 
specimen were measured and recorded on Mosely recorders shown in figure 23. 
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HYDROGEN - FREON 21 CRYOGENIC HEAT EXCHANGER 
FIGURE 16 
fid 




FIGURE 18 
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Test Setup 


FIGURE 19 
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Test Specimen 


FIGURE 21 
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Strip Chart Recorder 
FIGURE 23 
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Freon Pump 

FIGURE 24 
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Hydrogen inlet temperature was sensed with a Rosemount platinum probe and all 
other temperatures were sensed by calibrated copper constantan thermocouples. 
Hydrogen gas flow was measured using a Daniel's orifice as shown in figure 22. 
Hydrogen pressure and flowwere controlled by manual valves upstream and down- 
stream of the specimen. Hydrogen gas was discharged to a stack. Figures 19 
and 21 show the specimen installed in the test rig. Note that liquid droplets 

on the specimen god adjacent ft&ing were caused by a rein Viewer preceding 

photography, and were not present during testing. 

Freon 21 was circulated in a closed loop pressurized with nitrogen in 
order to insure that the Freon remained liquid. The Freon puny and ullage 
tanks are shown in figure 24. Flow rate was controlled by throttling valves 
and a puny bypass valve and was measured using a turbine flowmeter and 
recorded on a 4 channel recorder. The specimen was supplied with Freon 21 in 
a counterflow direction through one or both hot and cold heat exchangers 
via a mixing valve. Positioning the valve permitted assumption of the various 
heat loads as well as allowing specimen inlet temperature transients. Inlet 
and outlet Freon pressures were measured on 200 psi gages. Inlet and outlet 
tenyeratures were sensed with calibrated copper constantan thermocouples and 
were recorded on a f-fosely recorder. Differential pressure across the Freon 
side of the specimen was sensed by a pressure transducer and also recorded 
on a Mosely recorder. A detailedlist of instrumentation is contained in 
Appendix B. 


Conduct of Test 


In conducting the test the Freon loop was pressurized with dry nitrogen 
to the top of the inlet ullage tank. The Freon puny was then started and 
Freon flow and specimen inlet tenyeratures were adjusted. Hydrogen flow 
was next initiated through the dummy specimen until a minimum inlet tenyerature 
was achieved , thereupon the entire "flow was diverted through the specimen. 

Steady state readings were taken when Freon and hydrogen flows and temperatures 
were stabilized. Note that inherent specimen instability caused oscillation 
of hydrogen conditions so average readings were taken. Testing was conducted 
in the following sequence as described in the test plan contained in Appendix A 

1. Steady state 5000 Btu/hr. 

2. Freon tenyerature ramp change from minimum to maximum. 

3. Freon tenyerature step change from minimum to maximum. 

4. Shutdown and startup of hydrogen at 5000 Btu/hr load. 

5. Steady state 20,000 Btu/hr. 
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6. Steady state 35,000 Btu/hr. 

7. Steady state 50,000 Btu/hr. 

8. Hydrogen flow shutoff at 50,000 Btu/hr. 

9. Hydrogen flow startup at 50,000 Btu/hr. 
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STEADY STATE PERFORMANCE 


Steady state performance tests were made on a tube - in - tube heat 
exchanger sized for the projected Shuttle application in order to verify 
the computerized analytical design procedures. In the tests Freon flow 
rate and temperatures were set at values approximating those expected 
in Shuttle and the flow rate of cryogenic H 2 at supercritical pressure was 
set to obtain nominal heat rejection rates of 5000,. 20,000, 35,000 and 
50,000 Btu/hr. A summary tabulation of the steady state data is presented 
in Table II. TWo or three data points were taken at each nominal heat 
load and appear consistent. 

At heat loads of 35,000 and 50,000 Btu/hr oscillations in H 2 flow 
and pressure on the order of ±10% of total values were experienced. A 
description of the oscillations, together with a discussion of their 
possible impact on the Shuttle system, is presented following the transient 
test results. Since an average value for fluctuating H 2 pressures, inlet 
temperature and flow had to be estimated to evaluate steady state conditions, 
this added somewhat to normal experimental uncertainties. 

The percent error in the heat balance obtained between Freon and H 2 
is listed for each data point in Table II. Excluding the 5000 Btu/hr 
points, the maximum percentage heat balance error was 7.7% at a 50,000 
Btu/hr point, and consequently within acceptable experimental error. 

The three 5000 Btu/hr points have considerably higher percentage error, 
but because the Freon inlet to outlet temperature difference was small at 
this heat load an absolute error less than 2°R in a Freon temperature 
reading would account for the error at two of the data points . It appears 
there may have been a systematic heat balance shift between different heat 
loads, illustrated in figure 25, which has not been explained. Heat loss 
or gain to the Freon through the outer heat exchanger tube wall would pro- 
duce a shift in the indicated direction but is calculated to be relatively 
small. 

The values for flow rates , pressures and inlet temperatures , recorded 
during testing for one data point at each nominal heat load were input to 
the computerized heat exchanger model. The resultant predictions at these 
conditions are listed in Table II. The predicted pressure drops are for 
the heat exchanger tubes only. Consequently the experimental pressure drops 
which include fittings are higher. The experimental Freon AP was erratic 
throughout the testing and is not believed to be precise. 

The agreement between experimental and analytical values of the primary 
variable to be controlled by the heat exchanger, the Freon outlet temperature, 
is seen to be excellent. The maximum difference is 2. l d R at the 5000 Btu/hr 
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point, which represents a 14.7% difference between the relatively small 
experimental and predicted Freom AT in the heat exchanger. At the three 
other heat loads the differences between the experimental and analytical 
Freon AT are very small percentage wise. The experimental heat imbalance 
appears primarily as a difference between experimental and predicted H 2 
outlet temperature , suggesting the imbalance results from H 2 heat capacity 
or temperature readings. These could be the result of hydrogen oscillations 
at the higher heat loads. 


TRANSIENT TEST RESULTS 


The cryogenic heat exchanger would be employed for heat rejection from 
the Shuttle Freon loop in such a way that a change in the system heat load 
would result in a change of the Freon temperature entering the heat exchanger. 
The heat exchanger control must then respond to the changed heat load by 
changing H 2 flow to the heat exchanger so as to maintain the Freon exit 
temperature within a relatively narrow control band. 

In order to design an appropriate control the "open loop" response of 
the Freon exit temperature to a change in both Freon inlet temperature and 
H 9 flow rate, as determined by the thermal dynamics of the heat exchanger 
itself, must be predicted. An analytical model of the heat exchanger 
dynamics was programmed for a computer for the purpose of control analysis. 
This model was described in NASA CR 115569 and is included in the dynamic 
model of the Shuttle Freon loop described elsewhere in this report. The 
primary purpose of the transient testing reported below was to verify the 
accuracy of this analytical model. 

Tests were perfoimed in which there was: 

1. An independent change in H 2 flow into the heat exchanger (with Freon 
inlet temperature and other variables held constant^-*-)) , and 

2. An independent change in Freon temperature at the inlet to the heat 
exchanger. 


(1) Due to practical limitations in the test setup, it was not possible 
always to maintain exactly constant values of variables other than the 
one deliberately changed. For instance, a variation in Freon flow resulted 
from changing the Freon inlet temperature because the mixing valve for 
Freon from two different temperature reservoirs was not linear. However, 
such undesired changes were small enough that they were adequately com- 
pensated for in correlating the data by using average values over the 
period of the transient. 
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The resultant change in Freon exit temperature as a function of time 
was recorded (along with other significant variables) and compared to 
predictions obtained from the computerized model. Although with a control 
system in operation a Freon inlet temperature change would generally be 
accompanied by a simultaneous change in H 2 flow rate, it is desirable from 
an analytic viewpoint to study each separately as was done in the testing 
because the time response of the Freon exit temperature is significantly 
different in the two cases. 


Freon Exit Temperature Response to H 2 Flow Changes 


Figure 26 shows the recorded experimental response to shutting off 
H 2 flow at a beginning nominal heat load of 50,000 Btu/hr. This data has 
been plotted in Figure 27. Figure 28 is a plot of Freon exit temperature 
response data when H 2 was first turned off, then turned on at a nominal 
heat load of 5000 Btu/hr. In the tests, H 2 flow to the heat exchanger 
was regulated by a hand valve and an off or on "step" change in H 2 flow 
took place over the finite time required to turn the valve. Inspection 
of the data indicates this time was on the order of 2 or 3 seconds . 

Prediction of these Freon exit temperature transients obtained from 
the computerized model are also plotted in figures 27 and 28. An instantaneous 
step change in H 2 flow was input to the computer, so that the analytical line 
leads the data by about two seconds and does not show the rounding at the 
beginning of the transient produced by closing the valve over a period of 
2 or 3 seconds. However, the slope of the analytical Freon exit temperature 
response very closely matches the data. This may be seen most clearly in 
figure 27 where the transient extends over a sufficient range for the 
slope to be apparent. 

The nature of the response of the Freon exit temperature to a step change 
in H 2 flow may be seen in these plots. A change in Freon exit temperature 
begins almost immediately after H 2 flow is changed, and exhibits a nearly 
constant slope over the greater portion of the transient. These two char- 
acteristics result from the fact that the time constant of the heat exchanger 
wall separating the H 2 and Freon is small, both absolutely and relative to 
the transport time of the Freon in the heat exchanger. 


Freon Exit Temperature Response to Freon Inlet Temperature Change 

Figure 29 is a reproduction of all the recorded variables during a 
transient initiated by a change in Freon inlet temperature . In figure 30 
the inlet temperature change and the outlet tenperature response have been 
plotted from the data. It was possible to input to the computerized model 
a Freon inlet tenperature change very close to that experimentally recorded. 
As a consequence, the predicted exit temperature response plotted in 
figure 30 closely matches the data throughout the entire transient (unlike 
the H 2 flow change transients where the initial H 2 flow change could not 
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be exactly input) . The slight divergence at the end of the transient shown 
in figure 20 is caused by a shift in Freon flow Tate. 

The difference between the Freon exit temperature response to a change 
in H 2 flow and a change in Freon inlet temperature may be seen by comparing 
figures 27 and 30. The delay in response shown in figure 30 is due in large 
part to the time it takes a particle of Freon to traverse the heat exchanger. 

DISCUSSION OF INSTABILITIES ENCOUNTERED DURING TESTING 


Flow and pressure instabilities have been widely documented when 
cryogenic and nan-cryogenic fluids are vaporized under forced convection, 
and when supercritical fluids are heated. Two -phase flow in a fluid being 
vaporized in a heat exchanger is always unstable to some degree (if only on 
a microscopic level) because interrelations between heat transfer and 
fluid momentum are typically such that local perturbations in the density 
of the vaporizing fluid occur. The situation when a supercritical fluid 
is heated is analogous because gradients in fluid properties, both across 
the boundary layer and in the direction of flow, are similar to those in 
a vaporizing liquid. Under certain circumstances flow and pressure varia- 
tions in the heat exchanger become oscillatory and large in magnitude. 
Hypotheses concerning the "feedback" mechanisms which produce these oscil- 
lations will be discussed briefly further on. Such oscillations in 
pressure and/or flow are detrimental if undesirable effects are produced 
in the system of which the heat exchanger is a part. 

During the testing covered by this report the instrumentation employed 
was not designed to completely define the oscillations. For instance, the 
phase relations between inlet and outlet pressure variations were not con- 
tinuously recorded and would be of value. However, the following general 
description of the oscillations may be given. 

The maximum pressure and flow oscillations encountered were on the 
order of +101 of total values and were at a frequency of about one hertz. 

■pie oscillations had a limit-cycle nature: that is, they did not increase 

indefinitely but remained at the same amplitude. The largest amplitude 
appeared at the highest heat load (ft? 50,000 Btu/hr), while at the lowest 
load (ft# 5000 Btu/hr) none were detectable on total pressure or flow gages 
although variations appeared on the H 2 heat exchanger AP transducer and the 
H 2 inlet thermocouples. Flow variations produced variations in the thermo- 
couple reading H 2 inlet temperature, but often none in the H 2 outlet thermo- 
couple and never any significant variations in thermocouples reading Freon 
temperatures. Variations in H 2 inlet temperature were nearly linear (i.e., 
sinusoidal) while those appearing on the AP transducer were not. Recorded 
oscillations in H 2 inlet temperature and AP may be seen in figure 29. 
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The effect of such oscillations as were encountered during testing may 
be reasoned to be non - de tr iment al to the Shuttle system in which it is pro- 
posed to use the cryogenic heat exchanger for the following reasons: 

1. The instability was not of sufficient severity to cause oscillations 
in the Freon outlet temperature . Consequently there is no possibility 
of unstable feedback through the heat exchanger control which senses 
and controls Freon outlet temperature by modulating H? flow; the 
controls would not "know the H 2 flow was oscillating." 

2. The experimental steady-state and transient thermal performance of 
the heat exchanger appeared not to be affected by the oscillating 
flow, and is analytically predictable. 

3. The Shuttle H 2 feed system proposed for the cryogenic heat exchanger 
consists of a large storage tank supplying to fuel cells. H 2 

is to be returned from the cryogenic heat exchanger to, the tank and 
the heat exchanger serves the supplementary function of pressurizing 
the tank by adding thermal energy. This large tank may be expected 
to be an efficient damper for pressure oscillations originating in 
the cryogenic heat exchanger, so no detrimental effects on tank 
pressure control are expected. A "hard" H? circulating punp, i.e. , 
one with steep pressure versus flow characteristics, would prevent 
pressure signal feedback between the tank inlet and outlet and 
therefore prevent a possible resonant condition within the tank. 

Although no detrimental system effects are expected from such oscillations 
as were encountered in the testing reported, further definition of the nature 
of the instabilities encountered is clearly required before such a system 
could be designed with confidence , for at least two reasons : 

1. The system itself has not been completely defined and is subject to 
change as Shuttle design progresses. 

2, Since the tests were designed only to obtain steady state and tran- 
sient thermal performance of the heat exchanger, the interface bet- 
ween the heat exchanger and its H 2 feed and exhaust systems were 
not designed to simulate an indefinite H? system. H 2 pressure drop 
across the inlet valve was much larger than in the contemplated 
system, while a valve having choked flow at the heat exchanger 
exit used in the tests would not be part of the system. (The 
effect of the large inlet pressure drop would be to decrease the 
magnitude of pressure fluctuations in the system while the choked 
orifice would be expected to increase their magnitude.) In the 
absence of an accepted general theory concerning the mechanisms 
underlying instability in supercritical heat exchangers, data can- 
not be extrapolated with confidence to other conditions. 


87 



SVHSER 6180 


Hamilton 

Standard 


U 

DIVISION or UNITED AIRCRAFT CQ*W»OF*TlO'+ 


P® 


r 


A summary o£ literature on experimental and analytical investigations 
into instabilities in vaporizing and supercritical heat exchangers is pre- 
sented below to provide background. The discussion is then concluded with 
recommendations for further investigations aimed at the Shuttle application. 

It should be noted at this point that even without a cryogenic heat exchanger 
rejecting heat from the Freon coolant loop to the external pressurization 
loop for a fuel cell H 2 tank, instability may be encountered in this external 
pressurization loop when heat is added by an electric heater for the purpose 
of tank pressure control. 

The common starting points for formulating an anlytical model of 
oscillations in vaporizing and supercritical heat exchangers is of course 
obvious : 

Write the partical differential equations for fluid continuity and 
momentum, and for energy transfer to the fluid in the heat exchanger. 

A problem arises in the region where the fluid is vaporizing or super- 
critical because density gradients and gradients in other fluid pro- 
perties occur in a direction perpendicular to flow due to heat addition. 
Seme of the known potentially de- stabilizing factors in two-phase flow 
with heat addition, such as slip between phases in the direction of 
flow and a non-linear dependence of heat transfer on the wall to fluid 
(bulk) temperature difference are difficult to investigate experimentally 
and exact theoretical relationships to other variables have not been 
established. 

Because of the complexities of the situation, some of which are described 
above, an exact solution to the differential equations is not possible: 
extreme simplifications are required to obtain a manageable solution. Depend- 
ing on which features are retained in a simplified model somewhat different 
pictures of the mechanisms underlying oscillations are obtained. For example, 
the oscillatory storage of fluid mass in the heat exchanger that must occur 
for flow to oscillate may be envisioned as: 

1. Due to variations in the width of a liquid or heavy supercritical 
fluid core riding on a vapor film created by heat addition at the 
wall (6) ; 

2. due to the piston like movement of such a core; or 

3. due to fluid density waves traversing the heat exchanger (4) (7). 

The three pictures need not be exclusive. 

A commonly employed simplification is to linearize the differential 
equations so that Laplace transforms may be taken and stability criteria 
directly applied. Freidly and Krishman (1) have recently applied Poincare* s 
method to obtain a model providing limit cycle behavior, which the oscillations 
usually exhibit and which a linearized analysis cannot predict. 
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Critical parameters separating stable from unstable conditions which 
result from different models include: 

1. Density ratio between the fluid at the heat exchanger inlet and 
outlet (4) ; 

2. the ratio of heat flux to mass flow rate in the vaporizing or 
supercritical region (6) ; and 

3. ' the ratio of fluid pressure drops in the system upstream and down- 

stream of the point at which the principal density change may be 
thought of as occurring in the heat exchanger (e.g. , the end of the 
vaporizing section) (2) . 

Mich more experimental work is required to determine which simplifications 
can be most fruitfully employed to construct a model which can be used as a 
tool for predicting instability in the Shuttle application. In the absence 
of such a tool, experimental data collected under conditions most similar to 
the one being considered should prove most reliable. The most extensive data 
on oscillations resulting from heating supercritical H? in tubes has been 
taken by Thurston (6) , who had success in correlating nis own data using the 
parameters: 

heat flux 

NgO = mass flow - rate x heat oiTvaporization 

specific volume change due to vaporization 
N SV “ specific" volume of"” dense - phase 

At the inception of oscillations, 

% = .0045 (Ngv)'* 75 


This correlation must be applied with caution, however, when conditions 
are different. One difference between Thurston's tests and the tests reported 
here is that Thurston's tube was electrically heated, giving a different 
response of wall temperature to a change in the H 2 heat transfer coefficient 
than would occur with a fluid being cooled on the other side of the tube. 

Further investigations, both experimental and analytical, into the nature 
of the mechanisms producing oscillations, such as were encountered in the tests 
reported, would clearly be of value as a design guide in cryogenic heat 
exchanger applications. The sort of design changes that, might be made to 
stabilize the heat exchanger under the conditions tested are suggested by 
previous work. For example, using Thurston's criterion, stable H? flow 
could be obtained by insulating the tube wall in the supercritical region 
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and having the H 2 flcsw passage provided by a spiral insert inside the tube 
in the supercritical region. However, such measures cannot be taken with 
confidence without better understanding. 

Such tests might be most productive if conducted with scaled-down heat 
exchangers. In spite of the obvious danger of scaling when the phenomenon is 
not completely understood and the difficulties in accurately instrumenting 
the tests in a miniature scale , the cost of employing the sophisticated 
instrumentation that could be required together with large H 2 flows would appear 
prohibitive. 
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1.0 GENERAL INFORMATION 

1.1 Scope 

The purpose of this plan is to define the procedures and equipment 
to be utilized to test a tube - in - tube cryogenic heat exchanger. 

The objective of the test program is the correlation of analytically 
determined steady state and dynamic values witd test results, Test- 
ing will consist of steady state runs for thermal performance and 
transient runs for evaluation of the transfer function. Performance 
will be evaluated at four different conditions and transient operat- 
ing characteristics will be measured at five conditions . Wyle 
Laboratories of Norco, California will conduct all testing utilizing 
super critical hydrogen and Freon 21. 

1.2 Applicable Docunents 

NASA/MSC Contract NAS 9-12725 

Hamilton Standard Program Operating Plan B88-001 

Hibe - in - tube heat exchanger drawing SVSK 86021 

1.3 Functional Requirements 

■pie tube - in - tube counterflow heat exchanger is designed to operate | A 
in the Space Shuttle Orbiter. Stored sub critical hydrogen provides 
a nominal 50,000 Btu/hr of cooling. The heat exchanger shall be 
capable of operation at heat loads as low as 5000 Btu/hr without 
freezing. 
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2.0 


2.1 


2.2 


2.3 


TEST PLAN AND PROCEDURES 


The heat exchanger will be subjected to steady state tests to 
determine thermal performance and transient tests for dynamic 
analysis. 

S teady State Tests 

These tests will be conducted at each of the following heat loads: 

1. 5,000 Btu/hr 

2. 20,000 Btu/hr 

3. 35,000 Btu/hr 

4. 50,000 Btu/hr 

Transient Tests 

The following transients will be imposed on the heat exchanger: 

1. Freon inlet temperature change from 
52°F - 160°F at a uniform rate of 
20°F + 5°F/min. 

2. Freon inlet temperature change from 
52°F - 160°F in a step change. 

3. Start hydrogen flow @ 5,000 Btu/hr. 

4. Start hydrogen flow @ 50,000 Btu/hr. 

5. Stop hydrogen flow @ 50,000 Btu/hr. 

Test Conditions 


Conditions applicable to the above test runs are as follows: 


FREON HYDROGEN 


HEAT LOAD 

FLOW 

INLET 

OUTLET 

FLOW 

INLET 

INLET 

RATE 

TEMP. 

PRESS. 

RATE 

TEMP. 

PRESS. 

(Btu/hr) 

(lb/hr) 

(°F) 

(psia) 

(lb/hr) 

C°F) 

(psia) 

5,000 

1600 

52 

110 Min 

3.0 

-420 

190 Min 

20,000 

1600 

88 

110 Min 

12.0 

-420 

190 Min 

35,000 

1600 

124 

110 Min 

21.0 

-420 

190 Min 

50,000 

1600 

160 

110 Min 

30.5 

-420 

190 Min 


A- 3 



Hamilton 

Standard 


U 

QtvtSOhi OF UMITED *IAC#»AFT COlWCAiTifh 

fl® 


SVHSER 6180 


2 . 4 Test Procedure 


Definition of procedure, instrumentation and equipment are described 
in Wyle test procedure #3474 contained in Appendix A. 

2.5 Data 

The test data acquired during the steady state runs will be 
compared with prior performance calculations or will be recal- 
culated if data is not coincident with specific data points. The 
data obtained from the transient runs will be input to the computer 
for determination of effect on accuracy and stability. 

2.6 Test Schedule 


Testing will be started within 2-4 weeks of approval and will 
last approximately 2 weeks. 
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• EL SEGUNOO. CALIF. 


PURPOSE 

The purpose of this document is to present the procedures 

TO BE EMPLOYED DURING THE PERFORMANCE OF A TEST PROGRAM 
WHICH WILL BE CONDUCTED TO DETERMINE THE TRANSIENT AND 
STEADY STATE HEAT TRANSFER CHARACTERISTICS OF ONE HEAT 
Exchanger, Part Number SVSK86021 , 

REFERENCES 

Hamilton Standard Purchase Order Number SS043657NL. 

Hamilton Standard Test Plan AE-2-38. 

Hamilton Standard Drawing Number SVSK0G021 , dated 7-20-72, 
TEST CONDITIONS AND EQUIPMENT 
Ambient Conditions 

Unless otherwise specified herein, all tests will be per- 
formed at an atmospheric pressure of 29.92 ± 0,5 INCHES 
OF MERCURY ABSOLUTE, A TEMPERATURE OF +70 ± !5°F, AND A 
RELATIVE HUMIDITY OF APPROXIMATELY '50 PERCENT. 

Test Media 

All tests will be performed using Supercritical Liquid Hydrogen 

(190 psia min) (Inner Tube) and liquid Freon 21 (Outer Tube) as 2 

the test media. Freon outlet pressure shall be a minimum of 110 psia. 

Instrumentation and Equipment 

Measuring and test equipment, utilized in the performance 
of this contract, will be calibrated by the Wyle Labora- 
tories Standards Laboratory, or a commercial facility, 
utilizing reference standards (or interim standards) whose 

CALIBRATIONS HAVE BEEN CERTIFIED AS BEING TRACEABLE TO THE 

National Bureau of Standards. All reference standards 

UTILIZED IN THE ABOVE CALIBRATION SYSTEM, WILL BE SUPPORTED 
BY CERTIFICATES, REPORTS, OR DATA SHE ET S A TT E S T I NG TO THE 
DATE, ACCURACY, AND. CONTROL REQUIRED BY THE SUBJECT CONTRACT. 
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TEST CONDITIONS AND EQUIPMENT (CONTINUED) 

Instrument at ion and Equipment (Continued) 

Wyle Laboratories attests that the commercial sources pro- 
viding CALIBRATION SERVICES ON THE ABOVE REFERENCED EQU I P- 
MENT, OTHER THAN THE NATIONAL BUREAU OF STANDARDS, ARE IN 
FACT CAPABLE OF PERFORMING THE REQUIRED SERVICES TO THE 
SATISFACTION OF THE WYLE LABORATORIES QUALITY CONTROL DE- 
PARTMENT. Certificates and reports of all calibrations 

PERFORMED WILL BE RETAINED IN THE WYLE LABORATORIES QUALITY 

Control files and will be ava i ladle, for inspection, upon 
REQUEST, BY CUSTOMER AND/OR NASA REPRESENTATIVES. 

The following test equipment, or equivalent, will be used 

DURING THE PERFORMANCE OF THIS PROGRAM. 


EOU I PHENT 

Pressure Gauge 

Pressure Gauge 

AP Transducers (2) 

Temperature, Flow & 
Pressure Recorders 

Flow Meters (2) 


Mfg. & Model No. 
Ashcroft, 1279D 
Ashcroft, 1279D 
Validyne, DPI 5 
Moseley, 7100B 

Fo"xbo ro 


Range 
200 p s I g 
500 p s i g 
As Required 
As Required 


Acc y. 
±0.5$ FS 
±0.5$ FS 

N . A . 

±0.5$ 


As Required ±0.5$ 


The temperature probes will have a time constant of less than 0.1S 
seconds in water. 

The turbine flow meters will have a time constant less than 2S milliseconds. 


1 


Recorder time constant shall be less than 0.5 seconds. 

System Inspection and Test Witness 


Authorized Hamilton Standard and NASA representatives will 

HAVE THE RIGHT TO INSPECT THE TEST SYSTEM FOR COMPLIANCE 

w,i th Reference 2.1 and to witness all tests to be performed 
IN ACCORDANCE WITH THIS PROCEDURE. 
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4,0 


4,1 


4,1 .1 


4,1,2 


4.1.3 


4.1 .4 


4.1 .4.1 


4. 1.4. 2 


ffiQCEQURES A NO REQUIREMENTS 
Heat TranspjlrTes t 
Reference 2,2 


The Heat Exchanger will be installed in the Flow System 
shown in Figure 1 and subjected to liquid hydrogen and liquid 
Freon flow and temperature conditions in accordance with the 
steady state steps outlined below. 


2 


Freon Hydrogen 


Step 

No. 

BTU/Hr 

Flow 
Rate 
(Lb/Hr ) 

Inlet 

Temp, 

m 

Outlet 

Temp, 

F) 

Flow 
Rate 
,( Lb/Hr ) 

Inlet 

Temp. 

(°f) 

1 

5,000 

1600 

' '+52' 

+40 

3.0 

-420 

2 

20,000 

1600 

+88 

+40 

12.0 

-420 

3 

35,000 

1600 

+125 

+40 

21.0 

-420 

4 

50,000 

1600 

+1 60 

+40 

"'30.5 V. 

-420 


1 

2 


The test described in Paragraph 4.1.1 will be repeated ex- 
cept THAT THE FREON INLET TEMPERATURE TRANSITION FROM +52 |1 

TO +160 °F ^ WILL BE ACCOMPLISHED AT A RATEOF 20 ± 5°F PER l 

MINUTE. 


The TEST DESCRIBED IN PARAGRAPH 4.1.1 WILL BE REPEATED EX- 
CEPT THAT THE FREON INLET TEMPERATURE TRANSITION FROM +52 | 1 

TO +160V WILL BE ACCOMPLISHED IN ONE STEP AT THE MAXIMUM 
SYSTEM CAPABILITY. - 


During the steady state runs described in Paragraph 4.1.1, 

FREON OUTLET TEMPERATURE TRANSIENTS WILL BE DETERMINED AS 

follows: 


After Step 1 steady state conditions have been established, 

HYDROGEN FLOW WILL BE STOPPED FOR APPROXIMATELY 5 MINUTES. 

The hydrogen flow rate of 3.0 pounds per hour will then be 
instantly established and maintained until the freon outlet 
temperature reaches +40 0 F, 

After Step 4 steady state conditions have been established, 

THE HYDROGEN FLOW WILL BE SHUTOFF FO R AP P RO X I M A T E L Y 5 

minutes. The hydrogen flow rate or(3d.5i pounds per hour l| 1 2 

WILL THEN BE INSTANTLY ESTABLISHED AND MAINTAINED UNTIL THE ’ 
FREON OUTLET TEMPERATURE REACHES +40°F. THE HYDROGEN FLOW 
WILL THEN BE STOPPED AND THE FREON OUTLET TEMPERATURE WILL 
BE RECORDED UNTIL STABILIZATION OCCURS. 
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